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SHOPS IS 


In the present investigation, an attempt is made to 
understand the basic nature of gas-part iculate flows under 
various situations. Various problems that generally arise in 
Aeronautics, Industrial and Bio-fluid mechanics and natural 
aerodynamics are considered. Mathematical models based upon 
Navier -Stokes equations or its truncated form viz. boundary 
layer equations are considered* Particulate phase being 
pseudo fluid lacks definite boundary condition on the surface, 

A discussion of the various boundary conditions is undertaken. 
Various analytical techniques and numerical methods are used 
to solve these equations. Finally the global stability and 
uniqueness criteria of an IBVP (initial boundary value problem) 
of a dusty gas model are discussed. 

In Chapter I, basic assumptions involved in a mathematical 
formulation of a gas -part iculate flow problem are stated and 
a critical assessment of the governing equations namely Havier— 
Stokes equations and the derived boundary layer equations is 
included. A description of various boundary conditions on 
the surface for particulate phase is also incorporated. 

In Chapter II, gas-particulate boundary layer equations 
valid for flows on a flat plate are integrated numerically with 
and without the use of compatibility conditions on the surface. 
It is found that the surface characteristics in either case 



remain the same hut the detailed structure of the flow in the 
boundary layer is predicted better when the compatibility 
conditions are dropped from the analysis. 

In Chapter III, the problem related to the gas-particulate 
flows in between porous parallel plates with moving upper 
surface is considered. Both steady and unsteady problems are 
dealt with. Closed form solutions using various mathematical- 
techniques are derived. This study enhances our understanding 
of the basic nature of the flow in certain simple situations. 

In Chapter IV, the gas— particulate flow in tubes of 
varying cross-section is considered. Here, an asymptotic series 
solution of Havier -Stokes equations is developed. Three terms 
of the series are evaluated analytically. Zeroth order term 
gives the Poiseuille flow. Higher order terms give correction 
due to varying cross-section. For a diverging channel, the 
particles move towards the wall. Shear stress ait the wall 
decreases from its value for gas phase when particles are 
present. The present analysis is essentially valid for low 
Reynolds numbers. In case, we compute the flow for high 
Reynolds number, eddies are formed near the wall. 

In Chapter V, the problem of gas-particulate flow through 
curved pipes is solved. Dean (Mathemat ika, 1959, Vol. 13, p .77) 
discussed the motion of the gas phase in the curved pipe and 
developed analytical solution by using the method of separation 
of variables. This analysis Is extended to gas-particulate 



flows and the method of Green 1 s function is used, to solve 
the governing equations. The problem is solved for pipes of 
rectangular and circular cross-sections. It- is found that the 
addition of fine dust particles decreases the volume flow rate 
while the presence of coarse dust particles increases the 
volume flow rate. 

In Chapter 71, the global stability and uniqueness of an 
IB7P of a dusty gas model are discussed. Two types of uniqueness 
theorems of an IB7P for Navier -Stokes equations are pointed 
out* First of them ensures only one solution for each initial 
data while second guarantees only one steady state solution 
for not too small a viscosity. 



CHAPTE R I 
INIROEUGIION 

1.1 Genera l Introduction 

An understanding of the basic phenomenon in gas~p articulate 
flows is essential in Natural and Industrial Aerodynamics, 
several branches of engineering, environmental, biological 
and physical sciences etc. In two phase flow problems, 
suspended matter may consist of solid particles, liquid droplets, 
gas bubbles or some combination of these. In environmental 
sciences, we are greatly concerned with the dispersion and fall 
out of pollutants in air and the problem may be related to the 
nuclear fall out. Some of the other areas when the problem 
under consideration may find its applications are : movement of 
dust laden air; fluidization; use of dust in gas cooling systems 
to enhance heat transfer processes; coalescence of small 
droplets to form raindrops which might be considered as solid 
particles for the purpose of examining their movement; nuclear 
reactors with gas solid feeding; ablation cooling; batch 
settling; powder technology; dust collectors; sedimentation; 
acoustics; aerosol and paint spraying; blood flow with red cell . 
particles; pneumatic transport; slurry pipe lines; erosion of 
material due to continuous impingement of suspended particles in 


air etc • 
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During the past two decades considerable progress has 
been made in the development of continuum theories of 
particulate suspension. In a ccsn.tinu.um model, the particulate 
phase is considered as a cloud of particles* This cloud of 
particles is called pseudo fluid. The concept of pseudo 
fluid enables us to deduce the governing equations for gas- 
particulate flows similar to IT avier -Stokes equations for pure 
gas phase. In the present thesis, gas -particulate flows, 
flows with suspension and two phase flows have the same 
connotation. 

In § 1.2, we derive the governing equations for the 
two phase flow considering the gas and particle phases as 
compressible. Here, the terminology of compressibility for 
the particulate phase has been used only to indicate that 
the particulate density is a variable i.e, it is a function 
of space and time coordinates. Soo (1967) has called it as 
special compressibility. Vfe also deduce here the corresponding 
energy equations for the two phases. In order to gain some 
basic understaa ding of the underlying physical phenomenon of 
the two phase flows, laminar gas-particulate flows are 
considered in the thesis. Hence, the governing equations 

are also derived for laminar flows . 

' 1 . ' . . . . ' ' ' ■ ' 

In § 1.3, we have deduced the boundary layer equations 
when the particle phase is compressible and there is negligible 
volume fraction of the particles. 


In § 1.4, a discussion of 
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the boundary conditions is included. Special attention is 
given to the derivation of compatibility conditions for 
particulate phase , its physical interpretation and its 
drawbacks. Various forms of the governing equations so 
deduced are used in solving various physical problems in 
the following chapters* In § 1.5, a brief description of 
the various problems solved in the thesis, the method of 
solution and the important results obtained is included. 


1.2 Disc u ssion of the_ g c rver •ning_ equations of_ JQ o t ion^ of 
two J p~ha se""f Tow's ” 


Several investigators have derived the basic equations 
of the two phase flows. Significant contributions in this 
direction had been made by Soo (196 7), Marble (1962), 

Pai (19 70), Murray (1965), Van Deemter and Van Der Laan 
(1961), Drew and Seigal (1971) and Hinze (1963). Marble 
(1962) used the kinetic theory approach to derive the 
governing equations of motion. Van Deemter and Van Der 
laan (I96I) derived the continuity, momentum and mechanical 
energy equations in a formal way. They did not discuss 
the stress tensor and interaction forces. Conceptually 
their equations add little to our understanding of the 
basic nature of the flow of suspensions. Murray (1965) 
discussed the equations for two phase flows by considering 
the contributions of shear stresses in both the phases 
separately. His volume element of the mixture was made up 
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of fluid and pax-tide volumes. Murray (1965) derived these 
equations for fluidized beds but these are valid for dispersed 
two phase flows also. According to Drew and Seigal (19 71), 
equations derived "by Murray ( 1965) are not consistent as the 
equations are not valid for an arbitrary volume element of 
the fluid. Drew and Seigal (1971) derived the equations 
for two phase flows by averaging the various flow quantities 
over a given volume. These derivations are quite general 
are valid for a large class of problems. In a later paper. 
Drew and Seigal (1971) applied the averaged equations to the 
fluidized beds and foams. Pai (19 70) also used the 
phenomenological approach to derive these equations. He did 
not consider the effect of volume fraction of the particles 
on the stress tensor and the pressure gradient term. 

Here, we discuss the equations given by Hinze (1963) 
for two phase flows. These equations contain the effect of 
volume fraction of the particles. Although in the thesis* 
we consider negligible volume fraction of the particles in 
different flow problems, but this procedure enables us to 
envisage clearly the various terms that may imbibe the effect 
of the presence of the particles. It gives a better feeling 
of the 'approximations that have been incorporated by taking 
negligible volume fraction of the particles in the flow. 

Also, it opens the possibility of further extension of the 
problems discussed in the thesis where the volume fraction of 
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the particles becomes an important feature of the physical 
problem . 

Consider the fluid with velocity u^, temperature 1 

o o 

and density p with a cloud of solid spherical particles 

o 

having the same radius 'a' * The particle cloud is also 

described by the corresponding continuum variables u , 1 

P P 

and Pp . Here p = nm where n is the number of particles 
per unit volume and m is the mass of each particle . We 
assume that the particles are sparsely distributed and non- 
interacting. Plow around each particle is uneffected by the 
presence of other particles in the flow. let a denotes 
the volume fraction of the particles in a unit volume, p 

sg 

and p the densities of the material of the gas and the 
sp 

particles respectively, then 

P g = (1- «) P sg . (1.1) 

p = a p (1.2) 

p sp K ' 

The density and average velocity of the mixture are 

defined as follows: 

fjj = (!-“) » S g + ■ a o sp (1.3) 

I'm = (l-a) u g + a ^ (Volume velocity) (1.4) 

P M U M = ( 1 ” a ) P S g u g + a P sp u p ( Mass velocity) (1.5) 

Equations of continuity and momentum for the fluid and particle 
phases are the following: 
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fa + 6% (P 8“gi> - 0 


( 1 . 6 ) 


m + air ( p p u pi) = 0 


(1.7) 


at ( p g u gi) + a^' <‘’g u gi u g3 > = 

J 


„ S J2, + . iL, ft. ) 
dx ± dx. ^ ID ' 


-(i-a ) Pgg g 6 i& + I 1 . 


(1.3) 


|-r ( P U . ) + ( P U . U . ) 

at v p pi' ax. v p p i p d 


a P sp s ( lt9 ) 


Here, Einstein summation convention is used. In eqs . (1.8) 
and (1.9), £ is the direction parallel to the gravitational 

4- To 

acceleration 1 g ’ and represents the i component of 

the iter action forces due to virtual mass, buoyancy, slcin- 
friction resistance due to particles and resistance due to 
pressure gradient in the fluid phase. 

.Du . Du . 

E . = B a p ( --Jkt — - s a p r . „ 

•> sg v Dt Dt 1 & sg °i* 


Dt. 


where 


- R. + a'M 
i a XjL 


(l.io) 


D 

dt 


eft' + u gk ax 


k 


(1.11) 


D 

Dt" 


P 


and 


T . 

3-D 


. a , „ s 

M mmmsn# -+* 1 j msam.- mmxt 

5t pk 

v (Ik + . 8 hU) 

m v 0x. 6x i ; 


2 u 6U Mk . 

3 m 6 x^" °ij 


( 1 . 12 ) 


(1.13) 


u m being the viscosity of the mixture. 
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In eq. (1.8), Hinze (1963) lias introduced the stress 
tensor for the gas~particulate flow which incorporates the 
effect of the volume fraction of the particles in the 
description of the viscosity and the gradients of the 
mixture volume flow velocity. hue to the non interacting nature 
of the particles, stress tensor of the particles is not present 
in the momentum equation (1.9) for the particle phase hut due 
to the velocity of the particle phase, the mixture velocity 
is in general different from the velocity of either phase. 

Hence the deformation of the fluid element depends on the 
gradients of the mixture velocity. As is discussed earlier 
in the text, investigations of Murray (1965), Drew and Seigal 
(19 71) show that the stress tensor can he expressed differently. 


In the literature, various expressions for the coefficient 
of viscosity of suspensions ^ had been given. Einstein 
(1906) gave the following expression for the viscosity of the 
gas-particulate mixture 

“m = u g W + 2 - B “ ) (i-14) 

where v is the viscosity of the gas phase. Happel (1959) 

.'.O''.' 

gave the following expression 

— = 1 + 5.5a (1.15) 

y g 

I his formula differs from (1.14) because Happel considered 
the interaction forces of the particles. Eqs . (1.14) and (1.15) 
show that the viscosity of the gas increases by adding dust 
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particles to it. Inis fact was also supported experimentally 
by Ward and Whitmore (1950) and Williams (1953). 

Murray (1965) pointed out that static pressure in the 
particle phase is zero due to non interacting nature of 
the particles. Also, for a, mixture of practical interest, 
the number density of particles is insignificant in comparison 
to the number density of the gas molecules. Therefore the 
contribution of the particles to the pressure of the mixture 
is negligible. Rudinger (196 5) derived the following relation 
expressing the gas phase pressure in. terms of the mass 
fraction of the particles 


-Ji. - (1 + -Jba ?i-g)~ 1 (1.16) 

y v sp 4xa 

where p is the pressure due to the fluid phase, p the 

O 

pressure due to the mixture, 4 1 the mass fraction of the particles 

Iq 

and 1 = 2.69 x 10 cm is the Loschmidt number. He found 

that in case, the contribution of the solid particles to the 
pressure of the mixture is to be less than one percent 
(Pg/P > 0-99), then 


d > 0.0192 



, d = 2a 


(1.17) 


which shows that in most cases of interest, the particles 
contribution to pressure of the mixture is negligible if the 
particle diameter is less than a hundreth of a micron. He 
further suggested that the particles contribution to the pressure 
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can be neglected even at extremes of the density ratios. 


lor reasons stated above, equation of state is taken as 
p = R p . T , R being gas constant (1,18) 


T he term 


'O 

-D 


a p 


sg 



Du 


D~c 


PH 



in the interaction term (l.lo) corresponds to the force due to 
added mass of the particles. The numerical factor B 
corresponds to the volume of the virtual mass of the fluid 
subjected to the acceleration of the particles relative to 
the ambient fluid. Bor a spherical particle in an unbounded 
region, B = 0.5. Here B may differ from this value because 
of the presence of the neighbouring particles. 


Hinze (1963) did not add this term in the momentum 
equation (1.8) for the fluid phase but according to Newton's 
third law of motion, this term must also be added in the 
momentum equation (1.8) for the fluid phase. Murray (1965) 
pointed out that the form of this term is a conjecture. 

A term of the form 


Bap 


sg 


( V 


) + (% - u g ) 

. s i Pit g k 




K 


u 


can also be used. 

Inspite of the fact that pressure is negligible due to 
the solid particles, the term (-a |^“) 


in the solid phase is 
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introduced [Hinze (1963) (J mainly due to the resistance it 
offers to the particle motion. There is considerable 
controversy about the nature and the form of the pressure 
gradient term in the two phase flow equations. Some of 
these points are resolved by Shell and Soo (1979) and Boure (19 79) . 

R^ is the drag force due to difference in velocity in 
fluid and particle phases. lor the discussion of drag forces 
we shall consider the particles to be spherical in shape and 
of radius ’ a f . lor small R g (of order unity) Proudman and 
pe sirs on (19 57) gave the following formula for the drag 
coefficient 



where the leading term gives the Stokes law of resistance 
and the second term gives the correction due to the Stokes 
law by Oseen . 

In the literature, a large number of resistance laws 
incorporating the effect of Mach number, solid particle 
temperature, turbulence and rarefaction effects are available . 
Recently, Henderson (1976) gave such an emperical law, for the 
subsonic and supersonic regimes. Ilis formula for subsonic 
regime is valid up to Mach number M < 1.17. 

Por the sake of simplicity and ease in analytical 
calculations, we have adopted the Stokes law of resistance 
in our analysis. While passing, we may mention that the 
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free stream Reynolds number is large but the Reynolds number 
based upon the particle diameter and the relative velocity 
of the particle and the fluid may be sma.ll. Hence, the Stokes 
law may represent the resistance to the fluid with reasonable 
accuracy. If there are n particles in a unit volume of the 
fluid, then the resistance to the flow is 


E i = n 6 * 


which can be rewritten as 


where 


E i = 


R 


p R(u .-u .) 
p v gr pi' 


18 g 
-«r- “ 
d^P 

sp 


( 1 . 20 ) 


(1.21) 


Other interaction forces may be the lift and a torque 
acting on the particle . These forces would cause the particle 
translation in the direction normal to the flow as we 3.1 as 
particle rotation. Saffman (1965) calculated the lift force 
on a single particle. Rubinow and Keller (1961) gave an 
expression for the torque anting on a spherical particle. 

3? or small R , Hamed and Tab alt off (19*75) gave the following 
expressions for the forces acting along and perpendicular to 
the flow and torque which incorporated the results of Saffman 
(1965), Rubinow and Keller (I 96 I), Oseen and Tam (I 969 ) 


CCP 


R. 


^(u -u ) { 1+1 .5(2 a ) 1 / 2 + 3.75a } 

* m D 

2 ”|l/2 > 


x 


‘rji g P 
d 


{i+ f 5 r (u -u r+ ( v v r ~] 

16 V L v g P ' v g P ' * 4 

9 .69 ay„ 1 5 v„ Qu 


ltd 


(v -v ) P-- ( - «°~£) 

v g p ; L v ^ 5x ay ' 


( 1 . 22 ) 



12 


R 


a p 

— ? 


( v v ) {l+l ,5(2 a ) X / 2 + 3.75 a > 
rp S P 

{i + i [~;(u g -u p ) 2 + ( Vg ^ Vp ) 2 j 1/2 > 


9.69 a v r .,i dv dUp 1/2 

- s (V U P 5 U A 


T = 


7td 

.SR- 


a p _ 2 -i 5v„ 9u 


r 


Y j L>2 (5# “ -/') “ « 3 


(1.23) 

(1.24) 


In eqs. (1.22) to (1.24), (u g , v g ) and (u p , v p ) 

respectively denote the components of velocities of gas 

and particle phases for the two dimensional flow in x and 

y directions, u the angular velocity of the particles, y • 

J 

the radius of gyration of a spherical solid particle about an 
axis passing through its center, v the kinematic viscosity 
of the fluid and 


d 2p 


1 


- 


t .£ss> 

r 60 y« 


Here, and x are called relaxation times for 

particle translation and rotation respectively. 

Using this formula for resistance, they solved the gas 
particulate flow due to the impulsive motion of an infinite 
plate and found that a demixed region with no particles 
present developed near the plate . This is due to the 
particles migrating away from the wall. 
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get Pgg^-yj, I s the 'buoyancy force acting on the particles. 
Energy e qua tions 

Hinze derived, the equations conserving the mechanical 
energy to determine pressure. He did not derive the equation 
conserving thermal energy and the equation of state. Probably, 
for these reasons , Murray (1965) has stated that Hinze (1963) 
did not consider gas phase as compressible. Murray (1965) 
derived the equation representing the conservation of 
mechanical and thermal energy of the mixture and then using 
the momentum equations, got the energy equations for separate 
phases. Drew and Seigal (1963) deduced the thermal energy 
equation consistent v/ith their approach of using average 
quantities for the various variables involved in the problem. 
Pai (1970) got the energy equations for separate phases but 
did not consider certain interaction forces. In the present 
analysis, we have derived the energy equations f or separate 
species, taking care of the interaction forces as represented 
by virtual mass, buoyancy and pressure gradient besides the 
resistance force. 'this derivation is consistent with our 
approach adopted in deriving the momentum equations. We 
have also put these equat ions in terms of specific enthalpy 
for the gas phase and specific energy for the particulate 
phase. These forms are for ready use in solving specific 
problem involving compressible gas phase. 
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The corresponding energy equations for the gas and 
particle phases can he written as 


D r- /_ . 1 ..2 


Dt 


L Prr( e ^+ 9 u ^i) 


g v g 2 ~gi 
3 sg 


r~™ [ u. -(~ £ . . + t ■ • ) + q . 1 

8x . gi'* °13 i] ; ^-gD -” 1 

J 


( 1- a ) p g 6 . „ u. . + u . F . + Q 

V j b u x? gx pi 1 I 


SD 
(1.25) 


35F' t.Pp( e p + 2 U pi^ J u p i ap sp s6 i* ^pi^i - ^ C 1 * 26 ) 

Here, e and e denote respectively the specific 
£> P 

energies for g on and particle phases, q • represents the 

*o 3 

heat flux across the fluid element and is given hy 


8T s 

4 gj = - * 53“ 


1c being the coefficient of thermal conductivity of the fluid . 

In eq. (1.25), is the heat transferred from the 

particle cloud to the gas phase. l'he particles may have 
different temperature from the surrounding fluid. Due to 
this temperature difference, there will be heat transfer 
between the two phases. For the Stokes flow, heat transferred 
from the particles to the gas is given by 

p 0 (1 -I ) 

0 = 

IV 


where 


mO 

47tak 


is the relaxation time for temperature, G g is the specific heat 
of solid particles and 0 is the specific heat of fluid at 
constant pressure . 



15 


Multiply the momentum equation (1.8) by u • and then 
subtract the resulting equation from eq. (1.25). Using 
the continuity equation (1.6), we get the equation of 
energy for the fluid phase as 


D 

p g ift' 


5U 


xi 


+ (Qp+ 

J ci 


ID 6 


+ ( u P i- u gi) U 


= x = ~ X 

S PS P c 


(1.27) 

(1.28) 


I being the specific enthalpy of the fluid phase. Assuming 


G 


for the fluid phase; as constant, eq. (1.28) after 


differentiation and then using continuity equation (1.6) gives 


be 
p g b 


ST = o- 0. 


bf 


bp 


g P Dt ff 


Dt 


ir ~ p 


ax, 


(1.29) 


g g 1 

with eq. (1.29), eq. (1.27) becomes 

9 u o-i , 6 c -gi N 

p~ X rr- = + (X + 


J g p bt 


id ' 7? ax^ 

+ <v~ u gi> 1 


(1.30) 


Similarly multiply equation (1.9) by u. . and subtract the 

resulting equation from (1.26), using the continuity equation 

(1.7), we get the equation of energy for the particle phase as 
be 

P p E# = - % (1.51) 

Also 


C s *P 


(1.32) 
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with eq. (1.32), eq. (1.31) becomes 

DT 

P p °s IT# = “ % 


(1.33) 


Eqs. (1.6) to (1.9), (1.18), (1.30) and (1.33) form a 


set of seven equations for seven unknowns viz. 


u. 


'gi s ~p i 


u 


t a , 


p ,, p, 'i n . and T . For a particular flow problem, the above 
equations are solved with suitable boundary conditions. 

With rotating particles in the fluid, we need consider 
an additional equation representing the conservation of 
angular momentum of the particulate phase. Earned and 
fabakoff (19 74) gave the following equation for the conservation 
of angular momentum of the solid particles. 

n _ nw Y 2 

P P 3JT - T / Y j 

where ’l 1 is the torque of bulk interaction per unit volume 
between the particle and the fluid. 

Hamed and fabakoff (1975) further found that the presence 
of the solid particles with rotation causes an additional 
anti-symmetric stress tensor and gave the following expressions: 


xx 


2 

-p- g- 


/du 4. Hi J. p„ 5, u 
y g^ex + ay' 2y g 6x 


a = y (H + H) + a 

xy g K dy ax' 


2 i ai 

2 P p Y j EET. 


— H & ✓ du . o V\ . o 

- “ p_ 5 "gfe + 5y> + 2u 


yy 


av> 

yy‘ 


au 
'g ay 


a = y (hii + &JL) 
yx . g^ey dx' 


6U L ev>, a v 2 Du 

p Y. 


3 Dir 

J p 


where a a a and a are the elements of the fluid 
xx* xy* yx yy 
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stress tensor, (u,v) the components of velocity of the fluid 
phase in a two dimensional suspension flow. 


P artic ular f omn ._of_ the e _g ff _atiogas y>f motion 

We simplify eqs. (1.6) to (1.9), (1.18), (1.30) and 
(1.33) Toy making the following assumptions: 


(i) 


(ii) 


(iii) 

( iv) 


( v) 


Gas phase is incompressible i.e. p is constant 

o 

but the particle phase is compressible. 

The temperature difference bet ween the fluid and the 
particles is negligible. 

Volume fraction a of the particles is negligible. 
Body forces, rotation of the particles and the lift 
forces are negligible. 

The particles are sparsely distributed and spherical 


in shape. 

With these assumptions eqs. (1.14) and (l,15) give 

U m = y g = y ( sa y) (1.38) 

and the viscous stress tensor T . . (eq, (1.13)) becomes 

1 J 

(1.37) 


T. . = y ( „ JLsl) „ 3 . y £ 

id fix. ; 3 6x, 5 ij 

4- Jx 


The equations of continuity and momentum for the two 
p hases become 


9 

ax’ 
a p„ 


(u gi ) = 0 


+ 


ax, 


( P u , ) 
p pl ; 


= 0 


(1.40) 

(1.41) 
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Du - 

5T 


p ^,,£Uz 


9p 

dx i ' 5 x^ v ID 


+ ^ (x..) -FPpCUgi-^i) 


Du . 

p P Btp = ^Pp^gi^pi) 


(1.42) 

(1.43) 


Energy equations (1.30) and (1.33) become redundant in 
vie?i/ of the assumptions (i) and (ii). 


Here, the unknowns are u . ,p , p and u Eor a 

gx p P x 

2-D flow, the number of unknowns is six and we have six 
equations to determine them. Hence, with prescribed boundary 
conditions, it should be possible to solve a gas-particulate 
flow problem . 


Further, we notice that with the above assumptions, ' the 
gas phase realizes the presence of the particles only through 
the resistance term* Similarly, the particle phase 
equations are linked with the gas phase equations through the 
resistance term. Eor a particular flow problem, we have to 
solve the above set of equations simultaneously. 

1.3 L am inar bjeimdary^ la yer o f gas p fartjl^late^ flow s 

The concept of viscous flow boundary layer can be extended 
to gas-particulate f Iqws . The boundary layers with particle 

suspension may be turbulent. In the present study, laminar 
boundary layer of gas particle suspension is studied with the 
purpose of developing a basic understanding of the interaction 
of gas-particle flow near the solid wall. 


19 


The study of the boundary layer of gas-particulate flow 
is important for the following reasons: 

( i) To find the effect of the presence of the particles on 
the structure of the boundary layer and its surface 
characteristics such as skin friction. 

( ii) To study the particle retardation and accumulation 
on the surface. 


Here, we shall derive the boundary layer equations of the 
gas-particulate flow by on order of magnitude analysis for the 
incompressible gas phase and compressible particle phase, 
lor a 2-1' steady flow, the basic governing equations (1.40) 
to (1.43) are 


8u j. 9v - n 

dx + ay “ 0 


(1 .44) 


P ( u + v^) 

^ an ay' 


ap , „ /3 2 u a 2 u\ 


ax^ ay 


1 P (u-u ) 

p p 


P(u & + V ^V) 

^ ax ay y 


a v> 


2 2 

ap , ,, / a v , a vx 

ay + u ( :*T + “T' 
J ax ay 


(1.45) 


Ip ( v-v ) 
P^ p' 


a_ 

5x 


(p u ) +f—(pv) = 0 

' VA V\ / Pi TT V ' y\ n ' W 


P P 


ay v P p 


(1.46) 

(1.47) 


a u 


au 


P,( U n A#) = Fp (U~U ) 


V p ax 


p ay 


p- 


(1.48) 


p p ( u p ii, 

p p 6x 


a v 

+ V P ayb 


1 P ( V-V ) 

P^ P' 


(1.49) 
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Here the suffix ' g ! has been dropped from the variables 
for the gas phase. 

Based on the Prandtle's assumption, the effect of 

viscosity near a solid wall is realized only up to a small 

distance '6'. 1 6 1 is small in comparison to other 

dimensions involved in the flow problem. If x ^0(1), 

then y ^ 0(6) and the continuity equation (1.44) gives 

that v v o(S)» Prom the momentum equation (1.45), we 

2 

derive that y must be of 0(6 ) in order to retain the 
same order of inertia and viscous terms. Prom eq. (1.46) 
it can be easily shown that (^-v) is of the same order as the 
inertia terms. Since for the gas flow, inertia terms are 
of order 6, it gives ( v -v) = 0(6) • Hence o(v ) = 0(v ) 

ir ir 

= 0(6). Thus the boundary layer equations for fluid phase are 


du 

dx 


+ 


^ = 0 

07 


Pu 


d_u 

dx 


du _ 9p 

5 y g 



o 


_ M 
ay 


Bp (u-u ) 
P 


(1.50) 

(1.51) 

(1.52) 


Pol lowing the same reasoning, the corresponding equations 
for the particulate phase are the following: 


d_ 

dx 


(p u ) + 

v p p ' 


9_ 

ay 


( p V ) 

v p p ; 


= 0 


(1 .53) 


P„ u 
P P 


6U 

dx 


+ 


p v 
P P 


. u 

<L£. 

ay 


Pp p (u-u p ) 


(1 .54) 
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Q v 


p u 
P P 


E. 


fix 


■+* 


p V 

p p 


Lr 

ey 


F 


p (v- 

P K 


•V ) 

p 


(1.55) 


The eqs. (1.53) -to (1.55) are the same in form as 
equations (1.47) to (1.49). ’with "boundary layer concept, 
eq. (1,55) is of the order *6* while eqs. (1.53) and (1.54) 
are of order unity. With compressible particle phase 
assumptions, the unknowns for the particle phase are three 
viz. P p ,u p aJlc ^ v p * Hence, we need three equations to 
determine them. Thus, for mathematical consistency, we 
have to retain the equation (1.55) in the system of equations 
to describe the boundary layer flow. It is presumed here 
that retaining of lower order equation viz. eq. (1.55) will 
not effect the over all accuracy of the system of equations. 
It may probably improve upon it. 

Marble (1962), Soo (1967), Tabakoff and Earned (19 72) 
neglected eq* (1.55) in 'their formulation of the problem. 

They assumed v^ = v, then effectively reducing the number 
of unknowns by unity. Soo (I 968 ), Singleten (1966) and 
Otterman (1969) retained the equation (1.55) but kept v^ 7 ^ v. 

1 .4 Boundary condit ions 

Bor gas-particulate boundary layer flows, the boundary 
conditions for the gas phase arc the followings 

At y = 0,u=v=0 (1.56) 

As y -* u = u^ (1.57) 
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where u m is the free stream velocity at the edge of the 
boundary layer . 

for particle phase, the boundary conditions are: 

At y = 0, V ■ = 0 (1.58) 

P 

As y -*■ u = u and p = (1.59) 

j j p poo p poo v ' 

where u ^ and p ^ are the velocity and density of the 
P P 

pseudo fluid of particles at the edge of the boundary layer. 


The particle phase being pseudo fluid lacks definite 
boundary condition on u^ on the surface. This condition 
is provided by evaluating the continuity and momentum 
equations for the particle phase at the surface. The 
analytical conditions so obtained are known as compatibility 
conditions. Bq. (1.54), when evaluated at the wall gives 


° u p 

-J 


~ , ^-—Foru =u - Fx 

Sx 'y=0 pw p°° 

The continuity equation (1.5o) at the wall gives 


IZ\ 

3 y 1 y=o 


0 (u = 0 at y = 0 for all x) 


3 v 


S in ce v^ = 0 ( v) , we can t ake y ^ 

eq, (1.53), evaluated at the wall gives 

p — - 1 


0 


( p u ) 

v n n J 


0 or -£2 = 

p°° 


ax p pi'y=o P noo lm 


u 


(1.60) 


(1.61) 


Conditions (1.6o) and (1.61) are termed as compatibility 
c ond it ions f or t he p art icle p hase . 

Compatibility conditions for the gas phase only were first 
used by Pohlhauson (1921) while solving Barman integral of 


■boundary layer equations. Over a period of time, the uso 
of these conditions in integral methods have given results 
that are resonabl'y in accord with the exact solution of the 
boundary layer equations. 

Higher order compatibility conditions can also be 
obtained by differentiating momentum equation and evaluating 
the resulting equation at the wall jSchlichting (1962) 3* 

Jain and Bhatnagar (1962) used these higher order compatibility 
conditions for boundary layer flows with pressure gradient 
and got reasonably accurate results. 

Researchers justify the use of these conditions ' in 
solving boundary layer equations by the various arguments,. 

A few of them are listed below: 

( i) Since in the boundary layer flow, we are concerned 
greatly with the condition near the wall, more boundary 
conditions should be used near the wall rather than at the 
edge of the boundary layer. 

(ii) Boundary layer equations are valid at every point 
inside the boundary layer. Hence, these equations when 
evaluated at the wall should provide an additional boundary 
condition . 

In gas -particulate flow problem, Soo (1967), Tabakoff 
and Home a (1972) used the compatibility conditions (1.6o) 
and (1.61) for the particle phase in the development of the 
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integral method. Marble (1962) in developing the series 
solution did not use the compatibility conditions (1.6o) 
and (1,61). The results of these two approaches to solve 
the same problem differ significantly both qualitatively 
and quantitatively. lor example, from Soo (1967) and 
labalcoff's and Homed analysis (1972), we find that surface 
particulate velocity varies linearly with the longitudinal 
distance x* = (^) and at x*=l, u nw = 0 and p 


pw 


pw 


tends 


to infinity. Marble’s analysis gives u^ w = 0 for all x 

up to the first order of the series solution. This 

differing nature of the results gave rise to the question 

of validity of the compatibility conditions for the particle 

phase. Also, one finds that imposition of the condition 

u = 1-x* in the analysis forces the solution to give the 
pw ^ 

result that u “ 1-x*. In view of these observations, it 
pw 

was found necessary to probe further the use of the 
compatibility conditions in the analysis of the boundary 
layer flow. 


1 .5 Au tho r 1 s_ contributi on 

In chapter II, gas-particulate boundary layer equations 
(1.50) to (1.55) are integrated numerically with and without 
the use of computability conditions. We use the finite 
difference scheme due to Granfc-Nicholson to integrate the 
governing equations (1.5o) to (1.54) with compatibility 
conditions and equations (1.5o) to (1.55) without compatibility 
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conditions. It is interesting to note that- loth the analyses 
give linear variation of the surface particulate velocity. 

The analysis without compatibility conditions predicts better 
structure of the flow. 

In chapter III, with a view to understand some basic 

characteristics of the gas-particulate flow, we obtain some 

exact solutions of the eqs . (1.4=0) to (1.43) in cartesian 
coordinates. Problems considered are: 

( i) the steady flow between infinite porous parallel plates 
when the lower plate is at rest and the upper plate moves 
with uniform velocity. The lower plate is subjected to a 
constant suction (injection) velocity while the upper plate 
is subjected to injection (suction) at the same rate 
(ii) unsteady flow between infinite porous parallel plates 
when the upper plate starts moving with the velocity increasing 
or decreasing exponentially with time. At time t = 0, we 
get the steady flow case as a particular case. It is observed 
that velocity of either phase increases as the suction velocity 
at the lower plate increases. In case, the upper plate moves 
with velocity increasing (decreasing) with time, skin friction 
at the lower plate increases (decreases; non linearly with time. 
It is further observed that when the upper plate moves with 
accelerating velocity, velocity profiles for the gas and 
particle phases are almost linear and the difference between 
the velocities of both the phases increases with time. In 
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case, the upper plate moves with velocity decreasing with 
time, the difference between the velocities of both the 
phases decreases with tine. 

In chapter IV, we solve Navier-Stokes equations (1.4o) 
to (1.43) with Incompressible fluid phase and compressible 
particle phase in cylindrical polar coordinates. lie 
physical problem solved is the gas-particulate flow in tubes 
of varying cross-section. Blasius (I9I0), Tanner (1966), 
lee and Rung (l9 7o) and Manton (19 71) solved the same problem 
for fluid phase only. Kaimal (19 77) extended the work of 
Manton for gas-particulate flows. He considered both the 
gas phase and particle phase as incompressible , He found 
analytical expressions for the first three terms of the 
series solution. His solution indicates the format ion of 
eddies even at Reynolds number of unity. 

In the present investigation, we modify Kaimal’ s analysis 
to compressible particle phase. We develop an asymptotic 
series in terms of the small parameter e, which gives the 
rate of variation of the cross-section of the tube with the 
axial distance and calculate first three terms. It is 
obtained that particles move radially towards the wall which 
results in excessive accumulation of the particles there. 

We further observe that the shear stress at the wall decreases 
by adding the coarse dust particles. The present analysis is 
strictly valid for low Reynolds number. In case, we compute the 
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results for high Reynolds numbers, eddies aro formed at the 
wall and the separation of gas and particulate phases take place. 

In chapter V, we solve eqs. (1.40) to (1.43) in cylindrical 
polar coordinates and solve the problem of gas-particulate 
flow in the curved pipe. Here, we assume that radius of- the 
circle in which central line of the pipe is coiled is large 
in comparison to the radius of cross-section of the pipe. 

Also, the secondary motion for the gas phase is replaced by 
the uniform stream. Dean (1959) solved the problem for 
fluid phase in the curved pipe and found that the volume 
flow rate decreases as the secondary flow increases through 
the pipe. In the present investigation, secondary motion 
for the particulate phase is calculated from the governing 
equations. We use the method of Green’s function to solve 
the boundary value problem so obtained. It is observed that 
in case the particles are fine, the volume flow rate for the 
clean gas is more than the volume flow rate for gas-p articulate 
flow. In the case of coarse particles, gas-particulate 
flow predicts more volume flow rate than the one predicted 
by clean gas. In both the cases it is obtained that the 
region where the primary motion is greatest, shifts outward. 

In chapter VI, we discuss the stability and uniqueness 
criterion of an initial boundary value problem of a dusty gas 
model described by eqs. (1.4o) to (1.43). There are certain 
investigations in this direction. Crooke (19 76) examined a 
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uniqueness criterion for a dusty gas model and fandapat 
and Gupta (1976) obtained a universal stability theorem 
for a dusty gas model. In both the investigations, authors 
considered constant number density of the particle phase, 

Ian dap at and Gupta (19 76) did not obtain any uniqueness 
criteria.! while Grooke (1976) did not bother about universal 
stability theorem. 

In the present investigat ion, we improve the universal 
stability theorem of landapat and Gupta (19 76) and as a 
consequence, we establish a uniqueness criterion. At one 
point, we depart from the analysis of Dandapat and Gupta 
to obtain a much improved theorem on universal stability. 

It ensures a wider range of stability region. We further 
generalise the analysis to cover the case of non uniform 
number density of particle phase. We get only the uniqueness 
theorem in this case. 


CHAP TER II 

GAS -PANICULATE BOUNDARY LAYER PLOY OH A 
PLAT PLATE* 

2 .1 In tr o duct ion 

In the present chapter, the problem of boundary layer 
for gas-particulate flow has been solved with a view to under- 
stand the structure of the boundary layer flow and the use of 
the various boundary conditions. This chapter consists of 
two parts. In part I, we solve the gas-particulate boundary 
layer equations by finite difference scheme due to Crank- 
Hicholson using compatibility conditions for the particle 
phase. In part II, the gas-particulate boundary layer 
equations are integrated numerically with and without the 
use of compatibility conditions. 

Part I 

The solution of the governing equations of gas-particulate 
boundary layer flow has been obtained by the series method 
and by the approximate methods. The results by the two 
approaches differ in several essential aspects. Marble (1962) 
developed the series solution for an incompressible gas p hase 
and compressible particle phase and solved the first two terms 
of the series expansion. His analysis indicated that the 

* The" contents of this chapter have been published in 
Refs . 0:5, 46j . 
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p article velocity on the surface remained zero upto the 
first order term of the series solution. Soo (1967) solved 
the same problem by momentum integral method using linear 
profiles for the gas and for the particulate velocities. 

Dak ale off and Earned (19 72) modified the procedure by taking' 
a fourth degree profile for the gas phase and a similar 
profile for the particle phase. These analyses showed that 
the surface particulate velocity decreased linearly with 
the longitudinal distance ’ x’ « Also, particulate velocity 
in the boundary layer decreased near the surface and then 
increased monotonically to attain the free stream value at 
the edge of the boundary layer. For some values of the 
governing parameters, it became negative inside the boundary 
layer. These features of the gas-particulate boundary layer 
seemed implausible. Hence, a need was felt to find exact 
solution of these equations. 

We solve the same problem by integrating the governing 
equations of the two phase boundary layers numerically by a 
finite difference scheme due to Crank-Eicholson . We find 
tljat the present procedure predicts a better structure of the 
flow than is given by approximate method of Tabakoff and 
Earned (1972). Here, the particulate velocity remains positive 
and increases monotonically from the surface value to its free 
stream value. Also, the present results fail at a much later 
station on the flat plate than the results of Tabakoff and 
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Hamed (1972). Graphs are drawn to illustrate the effect 

of the changes in non-dimensional parameter I’ on the 

coefficient of shin-friction 0^ and the particulate 

velocity u 
,P 

Here, the particles are assumed to be spherical in 
shape. They are supposed to be sparsely distributed so 
that a particle is not submerged in the wake of any other 
particle. for the mathematical simplicity, Stokes drag 
law is used. 

2 .2 Gov erning equ ations and b oun dary conditions 


The governing equations for a steady 2-dimensional gas- 
particulate boundary layer flow on a flat plate are derived 
in chapter I £eqs, (1.50) to (1,55) J, We further assume 
that v p = v. As such, these equations become 


3u + 

6x 

MT . 0 

ay 


• 

(2.1) 

pu ^ + 
dx 

PV = 

ay 

2 

- 

ay 

P P (u-u ) 

p A p' 

(2.2) 

a 

ax 

( p u ) 
W P p' 

+ ( P 

6y v p 

V) ' = 0 

(2.3) 

5u 

6-U 




U + 

p 6x 

■■ TT* .njbwit 

iy 

= i'(u-u p 

) 

(2.4) 

Hepe, x,y 

are the 

coordinates measured along 

and 


perpendicular to the plate length, u and u p the longitudinal 
velocity components for the gas and particle phases respectively 
and v the transverse component of the velocity for both 
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the phases, p and o are the densities of the gas and 

P 

particle phases respectively, v the coefficient of viscosity 
for the gas phase and 

r 18 u 

!:• = --g. ~ 

d p 

sp 

d being the diameter of a spherical particle and p the 

sp 

material density of the particles. 


Bounda r y 
I lo r gas ph ase 


( i) 

At 

y = o ; 

(ii) 

As 

y - °°: 


(2.5) 

( 2 . 6 ) 


Since u = 0 at y = 0 for all x, eo . (2.1) when 
evaluated at the wall gives 

fyly=0 = 0 < 8 * 7 > 

Evaluating eq_ . (2.3) at the wall and then using 
©q_s . (2.5) and (2.7), we get 


(i) 


3x ^ P p u p^y=G 0 


71 


( 2 . 8 ) 




Mr 


Evaluating eg. (2.4) at the wall and then using 
conditions given in (2. 5),. we get -i 

'h • 1 . * " 


du 

H | = 

3x Jy=0 1 


( i;L ) ir=n = W X (2.9) 

Conditions given in eqs. (2.8) and (2.9) are termed as 
compatibility conditions for the particle phase. 
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(iii) As y - u p = and P p = P poo (2.10) 

Here, is the velocity of loth the phases in the 

free stream and p the density of the particles in the 

poo “ c 

free stream. 

We take the initiaJL profiles from the approximate 
results of fabakoff and I-Iamed (19 72) and Jain and Ghosh (19 79). 



+ v 


(2.15) 
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where 


Boundary 

f - :;p-- 
e f sp 

conditions 100001116: 

I Bor, gas phase 


(i) 

0 

I! 

>> 

u = 

( ii) 



As y -* °°, 

u = 


(2.16) 

(2.17) 


(i) 

<L. 

(Pn^)L =0 

(2.18) 


ax 


* * y=o 



'auL. 




( ii) 


j 

= -! P 

(2.19) 


_3 X 

- 

y=o 


( iii) 

As 

y -> 00 , u = P =1 

* p p 

(2.20) 


Init ial p rofi les 


I Initial profiles for u,v,u and corresponding 

to the approximate results of Tabakoff and Hamed (19 72) are 
the following: 


u = 2n - 2 n 3 

-4 A / <7r,2 7 3 ,4 *v 

-b n + g. ( n -3^ + 3 n - n ) 

(2.21) 
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(2.23) 
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where 


and 


ft and 


and 


J'xju p r z ) 4 

p = 1 + — r ~ r — (l-6n + 8n ~ 3n ) 


(2.24) 


= I 
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> n 


P & 


/ p' dy» 


!L e .i= 

v 


lx 

z = ^ 

are calculated from the following equations 
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The first term on L.H.S . of eq. (2.25) is different 
than that of the corresponding term in eq. (44) of Tabakoff 
and Hamed’s (1972) paper, because of miscalculation there. 

We integrate eqs . (2.25) and (2.26) by Taylor l s series 

method under the boundary conditions 6 = 6 0 = 0 at z = 0 

„ „ „ 'P 

and calculate the profiles for u, v, u and p at 

L P P 

x = 0 .1 from eqs. (2,21) to (2.24). 


II Initial profiles corresponding to the approximate 
solution of Jain and Ghosh (19 79) are 


u = l-(l-n) 3 


(2.27) 

t = (| " 2 - 2" S + | n 4 ) 

d6 

cS 

(2.28) 

3 p = l-(l-a 2 )(!-„) 3 


(2.29) 

P p = l-(l-a3)(l-P) 5 


(2.30) 


where 

2 

p = ?■ , a c = u , a rr = P __ * A = ~Tr . 

2 pw ’ o pw* ^2 

u and being the -velocity and density of the particle 

phase at the wall. 


A, a 2 and a^ are calculated from the following 
equations 



(a 2 +.072)(a 2 -5.696)(a 2 +1.364)/(4a 2 +3) 2 = .A(|^ - |g a 3 ) 

CO 

(2.32) 


(2.33) 


a 2 a 3 = 2 “ £ ^ a 2" i " a 3 ^ 

We integrate eqs. (2.31) to (2.33) Toy Tayler's series 
expansion under the boundary conditions 

At x = 0> A = 0* ag = a^ = 1 
and calculate the profiles for u,v,u^ and at x = 0 .1 

2.3 Method of solution 


We integrate eqs. (2.12) to ( 2 .15) together with the 
boundary conditions (2.16) to (2.2o) by a finite difference 
scheme due to Gc ante— Nicholson . The various terms in eqs. 

(2.12) to (2.15) are replaced by the following differences: 


_ 1 — — 

u m+l/2,n ~ 2 ^ u m+l,n + u m 5 n^ 


(2.34) 

(2.35) 




ay 


_L « mm . m* 

m+l/2,n ~ ^ u m+l,n+l~ u m+l,n-l + u m,n+l~ u m,n-l^ 

(2.36) 


g2- ^ 

< 8 “=tW/2,n = ^7-yg (% + l,n + l- 2f Wl, 


n 


+ u , i _ -,+ u ^,--2u +u ^ i ) (2.37) 

m+l,n-l m,n J .-l m,n m,n-~l' v ' 

Here, the suffixes m,n represent the value of the 

4- y% 

variable at the grid point represented by m step in 

"t 

x-direction and n step in y-direction. 


Similar expressions for v, u , and their derivatives 
are used in eqs. (2.12) to (2.15). We difference eq. (2.12) 
at the point (m+ , n - •g) end eqs. (2.13) to (2.15) at the 
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p dint 


where 


( m + ‘|j n) t The resulting difference equations are: 


v 


la +l/2,n~ v m+l/2$n~l’"* ^~^ u m+l,n + 'V.+l^n—l 


^m,n~^n,n~l^ 


hVl in+ l+ BA 


-jn+1 r ^n^m+ljn 1 " Sa u m+ljn~l 


(2.38) 

(2.39) 


A — ~^ i/2 ,n > _ 
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4Ay 


2(»y) ? 


U 
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'm+l,n 


F / T 


_ _ . ±ait + + £. ( p 

21 *" (Ay ) 2 2 p m+l,-" p 


Ax 

0 = ~ a£ _ _ 1 ^ 

n 4 a y 2(Ay) 2 


+ P ) 

n p m,n 


h = 


1 - (S 2 . + u 2 ) _ -EiiZSjn 

2Ax m+1 > n m > n 4Ay 

^Ajn+l ~ u m,n~l^ 

+ ^ u m,n+l~ 2 u ni,n + \i,n-l^ 

~ | (p + ” ) («L n ~u ) 

4 p m+l,n p m,n m,n p m+l,n p m,n 

= D* (2.40) 


A* ^ 
^ % 


+ B* P 


+ 0* P 


m +l,n+l n p m+l,n n p m-+l,n~l 


A = 
n 


v 


m+l/_2.n 
4 Ay 


B* 

n 


u 


?m+l n 1 , - ~ v 

- ■ 'm+l/2,n~ v m+l/2,n-l ' 


Ax 


4Ay 


where 
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C* = -A* 
n n 


i)* = 

n 


u 

P 


Pr 


m,n ^m.n p m,n ,~ 


- (Y 

. — v rr 


m+l/2,n / - - 

„ “ (p n - P n 
4aY • F m,n+1 


. -n"<v r* 

+ XL u 


4 Ay”” ^ m+l/2,n" m+l/2,n~l) 

) 

** (2.41) 


A)!* iT . _ 

p m+l,n+l n Pm+l,n n p m-i-l,n-l n 


+ CL u 


= XL 


where 


^ = v - 


B 


-*•* 

"h 

j 

n 


m+l/2,n/( 4s y) 

^m+l ,n P_p_ 
4X 2 


‘C 


ifr* .. 1 / *-2 


i3 = 
n 


(u + a 

2 a 5 p m+l,n p m,n 


) 


- >+§&£ (u ) 

4&y p m,n+l p m,n-*l 

+ (% + i n + %n _U p ^ 

’ ’ L m,n 


Boundary conditions (2.16) to (2.20') become 

%,!.= \1 = 0 for m l 1 


u m H = = p n =1 for m > 1 

’ p m,N p m,ll . 


u = u ~ P p Ax for m > 1 

p m+l, 1 p m,l 


Pp 


u 


u 


m+ 1,1 p m+l,l p m,l p m,l 


for m > 1 


(2.42) 

(2.43) 

(2.44) 

(2.45) 
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where n = 1 represents the plate surface and n = 1 the 
outer edge of the houndary layer . 


With prescribed initial values [eqs. (2,21) to (2.24) and 

eqs. (2.27) to (2.3o)[] and boundary conditions (2.42) to 

(2. 45), eqs. (2.38) to (2.41) are solved recursively at each 

step in x-direction , The set of eqs. (2.39) for n = 2 to 

N-l is solved first by matrix inversion method. In the 

beginning, the values of u „_ l1 : u t p and 

m+i,n P m+1?n P m+ i jn 

^m+1/2 n occuri21 g 1x1 tlie coefficients A n , B n , 0 n and I> n 

at the (m+1) step are replaced by the corresponding values 
~b Ti 

at the m step. The set of eqs. (2.39) is thus rendered 

linear and is solved by matrix inversion method. In later 


iterations, u m+ j_ n in. the coefficients are replaced by the 
values obtained in the previous iteration, while the values 


of v_, /0 u and p in the coefficients 

m+1 /2.n> P m+ i, n P m+ i, n 

A , B , 0 n and D n are kept constant. Iterations are 
repeated till the difference between the values of u at 

Vi 

the grid points in (m+1) step in two successive iterations 
is less than lo” 4 . Generally three iterations are sufficient 


to get the converged values of u. Next the set of eqs. (2.38) 


for n = 2,...E for v is solved explicitly. 


The values of u, v thus obtained are substituted in 
the set of eqs. (2.41) for n = 2,... 1-1 and is solved implicitly 
by matrix inversion method. Similar iterative procedure 
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(used for u) is used to get the converged value of u. * 

P 

She set of eq_s. (2.40) for n = 2,... N~1 i s linear and 
is solved implicitly for p . This cycle of operation is 
repeated until the maximum difference in the successive 

—4 

values of a variable at each grid point is less than 10 . 

Similar operations are carried out at each step in x-direction. 


We then calculate the skin friction coefficient by the 
following formula 


G 


f 


tH 

2 0U ” 



.fH = fS+im 

W e € 


(2.46) 


As the integration is carried downstream, the step 
size ax is increased by .01 after every 25 steps till 
Ax = .1 is reached* Thereafter, it remains constant * 

Ay = .2 is kept constant. At some point in the downstream 
direction , the boundary layer edge lies outside the range of 
integration. At this point, the difference in the values 

, . mmti mm "bjd *4- Tq 

of u, u and p at h and (U-l) point becomes greater 
• P P 

—3 

than 10 and the range of integration in the y~direction is 
increased by lo steps (fig. 1). The integration is carried 
on an IBM 7044 computer of I.I.T. Kanpur. Depending upon 
the severity of parameters, it takes 3 minutes to 2o minutes 
of computer time . 


2.4 Discuss ion of the results 

We present here some of the characteristic features of 
gas-particulate flow in the boundary layer . The results of 
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the present analysis are compared with the results of Tabakoff 
and Hamed (19 72). lor comparison, following values of the 
parameters are used 

u co = 60 .96 m/sec . 

p = .9 752 leg ./m^ 

y = 1.5415 x lo" 5 kg./m 3 

P 0 „ = 80 1.0, 1602.0, 2403 .0, 8010.0 kg./m 3 
sp 

cc~ 0*1, 0.2, 0*3, 0«4 
d = 50 y , lOO y , 250 y 
1 = .3048 m 

In Pig. 2, we have compared the skin friction coefficient 

predicted Toy the present method with the results obtained by 

Tabakoff and Earned (19 72) . The present results are computed 

with the initial profiles given in eq_s. (2.21) to (2.24) and 

eqs. (2.27) to (2.3o) for p = 2403 kg./m 3 , a = 0.2 and 

sp 

d = lOO y . We note that both the initial profiles predict 
the same values of G^. We further notice that the velocity 
profiles for the two phases are also identical for the two 
different initial profiles used in the computation. We thus 
surmise that the change in the initial profiles has no effect 
on the flow structure downstream. In the following results, 
we have used the initial profiles provided by Jain and 
Ghosh (19 79). This figure also indicates that the present 
results predict a lower value of than the results of 

Tabakoff and Earned (19 72). Por the prescribed values of the 
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parameters, the present results are obtained upto x = IV. 1. c 
At x = 17.1, the surface particulate density becomes 
excessively large and the analysis fails at this point. 

For the same values of the parameters, the approximate 
results of Tabakoff and Hamed (19 72) hold upto x = 13,7. 

Beyond this point, oscillations in the boundary layer thickness 
and skin-friction are set in and the results fail thereafter. 
This shows that present results have greater range of validity. 


In Figs. 3 and 4, we have compared the profiles for 

u and Up from the present calculations with the corresponding 

results of Tabakoff and Hamed (19 72) for p = 2403 kg./m 3 , 

sp 

a = 0.2 and d = 100 y • Fig. 2 shows that at x = 6 .0, gas 

phase velocity profiles predicted by the present method and 

the approximate method agree reasonably well. Present method 

predicts velocity profile that lies slightly above the 

approximate method of Tabakoff and Hamed (1972). In figure 3, 

present results show that u increases monotonically from 

the value at the wall to its asymptotic value unity at the 

outer edge of the boundary layer while the approximate method 

of Tabakoff and Hamed (19 72) gives the velocity profile for 

u which first decreases from its surface value and then 
P 

increases to approach unity at the boundary layer edge. At 
x = 13.0, particle velocity profile becomes negative in the 


boundary layer. This variation of is not physically 

plausible. In this respect, present method seems to predict 
a correct u -profile. 
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In Figs. 5 and 6, we find "the effect of changes in 

R on 0^ and u^ profiles. Pig. 5 shows that 

increases as P increases. Pig. 6 shows that u. increases 

P 

as P decreases. 

Part II 

2.5 The problem of gas-particulate boundary layer flow is 

reconsidered to study the effect of the compatibility conditions 
for particle phase on structure and surface characteristics 
of the flow. References are available in the literature 
which use and do not use the compatibility conditions for the 
particle phase in the calculation of boundary layer flow. 

Marble (1962) developed the series solution for incompressible 
gas phase and compressible particle phase without the use of 
compatibility conditions. His analysis indicated that the 
particle velocity on the surface remained zero up to the 
first order of the series solution. Singleton (1965) extended 

the analysis with compressible gas phase and obtained series 

. * 

solutions valid near the leading edge and far from the leading 
edge of the flat plate. Soo (1968) obtained the series 
solution using compatibility conditions on the flat plate 
by assuming p^/p « 1. This assumption neglected the 
interaction term from the momentum equation of the gas phase 
and he could consider the Blasius solution for the gas phase. 

Rar from the leading edge where particulate velocity is zero, 
he incorporated the diffusion of the particles in order to 
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avoid "their accumulation on the surface. Ear from the 
leading edge, this analysis predicted that decreased 
in the boundary layer and attained the value less than the 
free stream value, Soo (1967) solved the problem for 
incompressible gas phase by momentum integral method using 
linear profiles for the gas and for the particle velocities, 
lab ale off and Harned (1972) modified the procedure by talcing 
the fourth degree profile for the gas phase and similar profile 
for the particle phase. Both of these investigators used the 
compatibility conditions for the particle phase and obtained 
a linear variation of the surface particulate velocity with 
the longitudinal distance ‘x* on the plate. Jain and Ghosh 
(1979) developed an integral method dropping the compatibility 
conditions and obtained the solution for all values of the 
longitudinal distance 'x* on the plate. The present authors 
attempted to analyse the effect of the use of compatibility 
conditions in solving boundary layer equations on the flow 
structure and surface characteristics. 

We solve gas -particulate boundary layer equations with 
and without the use of compatibility conditions by Crank- 
h'icholson scheme of finite differencing. We find that in 
both the cases, surface particulate velocity varies linearly 
with the longitudinal distance ' x’-. With compatibility 
conditions, particulate density inside the boundary layer 
decreases from its surface value to a value much below the 
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free stream condition and then increases to attain the free 
stream value while without the use of compatibility conditions, 
the density decreases monotonically from its surface value to 
the free stream value , In this way, the results without • j 
compatibility conditions are superior to the results with 1 / 
compatibility conditions. I'or other variables of the flow, 
the results are similar but differ slightly in magnitude. 

In the absence of experimental data, it is difficult to say 
which of the two results Is accurate, 

2 .6 Gov erni ng e quations amd_ boundary ^conjiitions 

Equations governing the 2-dimensional steady gas 
particulate boundary layer flow on a flat plate are the 
following 



8u 

+ Pj = 0 



(2.47) 
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= f 

P( U-U p ) 

(2.50) 
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- e ^n 




% 

JB, 

di 

+ v — E- 

p sy 

= f 

P (v-v p ) 

(2.51) 

Here 

V 

P 

is the normal 

component of velocity 

of the 


particle phase and the rest of the variables are the same as 
in (2 *11) , 
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For the analysis with compatibility conditions v = v 

is assumed and then eqs. (2.47) to (2.50) after replacing 

v by v are solved for u, v, and p^. Without 

compatibility conditions, we assume that Vp ^ v and 

eq_s. (2,47) to (2.51) are solved for u,v,Up,Vp and . 

For some unknown reasons, without using the compatibility 

conditions, we could not get the solution of eq_s. (2.47) to 

(2,51) when v = v. As stated in chapter I, § 1.4, it is 
P 

presumed that v = Vp = o(6) and as such its retaining and 
dropping will not change the results substantially. 

Boun d ary c onditions 

I Withn^ c ompa tibili ty condit ion s 


(i) 

At 

II 

o 

PI 

It 

<11 

II 

<11 

11 

o 

(2.52) 

(il) 

As 

y - oo u = u = p, =1 

J ’ P P 

(2.53) 

II With 

compat ibility_ condit ions 


(i) 

At 

*<1 

II 

o 

Ml 

PI 

II 

<11 

II 

o 

(2.54) 

(ii) 

‘”(p„uJ|_ =0 

(2.55) 


P \ -cr 

0 ux 

* * y=0 



fdu ”"] 



( Hi) j 

j 

= -F P 

( 2 .56 ) 


Jx J 

t y=o 


(iv) 

As 

y - u = Up = P p = 1 . 

(2.57) 


The initial profiles are taken from eq_s . (2.27) to 
( 2 .30 )at x = 0 *1 • 
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2 .7 


Solution of equations (2.47) to (2,51) without 
c ompat ib ilit y (TondYtx on sT~ * ----- - — 


v/e divide the range of integration into grid points as 
defined in § 2 .3 of part I and solve eqs. (2.47) to (2.51) with 
boundary conditions (2.52) to (2.53) by finite difference method 
due to Grank -Nicholson . The various terms are approximated 

as follows: 


u 


L m+l/2,n 2^ u m+l } n’’" u, m,,n) 
Qx m+l/2,n' 

~) m+ l/e,n= "y 


(.f)m+l/2.n ( u m+l,n~ u m,n^ 

Ax 7 ’ 

v 1 = ^(^+l,n+l~%+l,n-l +{i m,n+l“%,n-l 


‘‘ticm+l/Z.K+l/Z si (3 »+l,n+l/2“ U m,n + l/2) 


S y m+l/2,n+l/2 A y 


~i ^ u m+l/2,n+l~ u m+l/2,n^ 


(4l) . 

Sy m+l/2,n 


2(a y) 


2 ^ u m+l,n+l“ 2u m+l, 


n 


+ %+l,n-l +u m,n+l“ 2 %i,n + %,n-l^ 


(2.58) 

(2.59) 

) 

('2 .60) 

(2.61) 

(2.62) 

(2.63) 


Using (2.58) to (2.63) and similar e 2 cpressions for 

other variables v 3 u , v and p , eqs. (2.47) to (2.51) are 

^ X -jP 1 11 

differenced at the points (m+ g», n- ??)j (m+ g>,n); (m+ g-,n+ gr)l 

(m+ n+ -*) and (m+ n- -*) respectively. The resulting 

difference equations are: 
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T m+l/2,n v m+l/2 ,n-l~ u m+l,n-l~%,n~ a m J n-l^ 


\i u m+l,n+l' Sa u m+ljn~l 


(2.64) 

(2.65) 
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n Pm+l,n+l ^ Pm+l,n 


(2.66) 
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P* u 


n p m+l,n+l 


+ % u d ' = *£ 

^ p m+l,n n 


( 2 . 67 ) 


where 


u. + u. _ 
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8Ay p m,n+l p m,n p m+l,n+l p m,n+l 
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+ n+l + U m+l,n + u ta,n+l + %.,] 


^m,n+l -^m,n 


P** v + Q** y = R** 

11 p m+l,n ' n p m+l,n~l 


+ Q** y 
, TL P 


(2.68) 


where 


P** = — — (u + u + u + u ) 

n 8 ax p m+l,n p m+l,n-l p m,n p m,n-l 


y v „ _ 

+ p m+ l,n + p m,n + F p 

8Ay 4Ay . * 
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+ u + u 


8ax p m+l,n p m+l,n-l p m,n p m,n-l 


v v _ „ 

P itH -l , n~l _ + -Lil 


R n -/■ CU m+ l/2,n + v m+l/2,n-f " V P m;J1 _;f V P m;J1 ^ 


+ -~ C( 
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Boundary conditions in difference form Become the 
following: 

3 1 Without compatibility co nditions 



= v m,l = 

; = 0 
p m, 1 

for 

m > 1 

(2.69) 

= 

= u = 

p m,h 

p =1 

p m,N 

for 

m > 1 

(2.70) 


W ith c omp at ibility cond it ions 

\l = \1 = 0 for m * 1 


PU uL = p_ uL tor m > 1 
p m+l,l p m+l,l p m,l p m, 1 


(2.71) 

(2.72) 


Pm+1 l = ^n for m > 1 

m+± » x p m,l 


U , T = U = P : 

m > H Pm,H Pm,H 


1 for m > 1, , 


(2.73) 
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With prescribed initial values at Xj= .1, system of 

eqs. (2.64) to (2.68) with boundary conditions (2,69) to (2.70) 

is solved recursively at the point (x-vt ax). System of 

eqs. (2.65) for n = 2, , . .b-l for u is solved first. In 

the beginning, the values of u , 1 , u ; p and 

m+1 ’ n p m+l,n p m+l,n 

vL , / 9 „ occur ing in the coefficients A. B , C! , 1 at the 

Ul“r -L/ (£ jX.1 XI II XI II 

"til 

(m+1) step are replaced by the corresponding values at the 
"b ll 

m step. The set of eqs. (2,65) is thus rendered linear and 

is solved by matrix inversion method. In later iterations? 

values of n in the coefficients A n , B n , G n and ' I> n 

are replaced by the values obtained ±n the previous iteration 

while the values of v_,-i /0 . u and p in the 

m+i/^n P m+ i jn p m+l,n 

coefficients A^, B , C n and are kept constant. 

Iterations are repeated till the difference between the values 
of u at the grid points in (m+l) step in two successive 
iterations is less than lo”" 1 . Generally three iterations are 
sufficient to get the converged values of u. Next, we solve 
explicitly the set of eq_s. (2.64) for n = 2,.,.N for v. 
With values of u,v thus obtained, the set of eqs. (2.68) 
for n = 2,...U for v is solved next. The values of 
u,Vp thus obtained are substituted in the set of eqs. (2.67) 
for n = and is solved explicitly to get the values 

of u . In the end, we solve explicitly the set of eqs. 
(2.66) for- n = 1, . . .1-1 for . This cycle of operation is 
repeated till the values of the variables at different grid 
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points at (m+l)"^ step differ Toy less than 10~^ in successive 
cycles, and the calculation proceeds to the next step . At 
each step, skin-friction coefficient is calculated hy 

the formula given in eq. (2,46), 

A constant step size of ax = ,0l and. Ay = ,2 is used. 
Both the step sizes are kept constant throughout the integration. 
When the computation proceeds in the downstream direction and 
the difference in u, u^ and p at and (N-l)'*'* 1 grid 
points at any station on the plate exceeds a prescribed 

value, the range of integration in the y-direction is 
increased so as to get a smooth variation of the dependent 
variables. Initially 40 steps in y-direction are used and 
when the boundary layer edge is reached, this range is 
increased by 10 steps. Depending on the severity of the 
variables, computer time on IB.I 7o44 computer varies from 10 
to 45 minutes. 

Ihe solution of the equations with compatibility 
conditions is discussed in Part I of the chapter. 

2 .8 Pi s ous sion of the result s 

We compute the various flow variables for same values 
of the parameters considered in §2.4 of Part I and the results 
obtained from the analysis without compatibility conditions 
are compared with the results with compatibility conditions 
discussed in Part I . 



54 


In Fig. 7, we compare the results of the skin -friction 

coefficient predicted by the analyses with and without 

compatibility conditions with the results for clean gas flow. 

We find that the presence of particles in the gas increases 

the value of skin-friction coefficient. Also, predicted 

with compatibility conditions has lower values than 0^ 

without compatibility conditions. In Fig. 8 we compare 

values of u obtained from the two analyses. It is 

surprising to note that the results obtained without the use 

of compatibility conditions are exactly the same as the values 

obtained with the use of compatibility conditions. Both the 

approaches give a linear variation of u with the distance 
r D pw 

1 x' on the plate. At x = 17.1. u =0 and p becomes 

PW pw 

excessively large. Beyond this point, oscillations in 

"pp-profile are set up and the method fails to proceed further. 

In Fig. 9, we compare the values of p calculated from both 

the analyses. The graph illustrates that the values of p 

x pw 

calculated from both the analyses agree upto x = 13.0 and 

differ slightly thereafter. However, there is a sharp rise 

in the values of p" beyond x = 13.0. 

pw 

In Fig. lo, we compare u-profile by both the analyses. 
We find that u-profile obtained from the analysis without 
compatibility conditions lies slightly below the profile 
predicted by the analysis with compatibility conditions. 

11 illustrates that the u -profile obtained by using 

Jr 
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compatibility conditions lies slightly above the profile 
obtained without using the compatibility conditions. This 
graph further shows the entirely different behaviour of 
particulate density in the boundary layer* Relaxing the 
compatibility conditions, the analysis predicts better 
structure of p -profile', The values of p^ decreases 
monotonically from the surface value to the free stream 
value at the edge of the boundary layer . The analysis with 
compatibility conditions shows that p^ first decreases 
from the surface value and goes much below unity in the 
boundary layer and then it starts increasing and attains 
the value of unity at the boundary layer edge. Thus, the 
analysis without compatibility conditions predicts better 
structure of the flow than the corresponding results with 
compatibility conditions. 

In Rigs. 12 and 13, we give the variation of and 

with R when the compatibility conditions are not used. 
Rig. 12 shows that decreases as R increases and Rig. 13 

gives that p increases as R increases. 
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IG. 10. GAS PHASE VELOCITY DISTRIBUTION ACROSS 
BOUNDARY LAYER AT X=80 
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0 H AP 1E R III 

GAS -P ARII 0ULA1E FLOW' IN JElWM POROUS PARALLEL PLAINS 
3 . 1 In t r o due t i on 

Closed form solutions for unsteady gas-particulate flow 
problems had been obtained by several authors. For example, 

Liu (I966) studied the flow induced in an incompressible dusty 
gas by an infinite flat plate oscillating in its own plane. 

He obtained analytical expressions for the velocitjr profiles 
for the gas and for the particle phases, shear stresses at the 
plate and discussed mechanical energy dissipation. He indicate^ 
that thickness of the viscous diffusion layer was decreased by 
the presence of the particles and the shear stresses at the 
plate were higher than the shear stresses for the gas phase 
only. In a later paper, Liu (196 7) studied the flow induced by 
an infinite flat plate suddenly set into motion parallel to its 
own plane in an incompressible dusty gas. He gave approximate 
results for the velocity profiles in the form of series solutions 
valid for small times and for large times and calculated the 
corresponding skin-friction coefficients. Ilis analysis indicated 
that for small times, viscous diffusion layer grows parabolically 
like (vt) 1 / 2 ( v being the kinematic viscosity for the gas phase 
and t the time variable), independent of the presence of the 
particle phase. For large times, viscous diffusion layer also 
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grows parabolically lout as (vt)^/^ where v is the kinematic 
viscosity depending upon the viscosity of the gas hut total 
density of the combined gas phase and solid particles. 

Healy and Yang (_19 ) studied the oscillating two phase 

flows with viscous incompressible fluid and particle phases 

in a channel under the action of body and pressure forces. 

He showed that for the case of fixed body force acting on the 

fluid, the system tended to behave as a particle free fluid in 

equilibrium. It was also shown that the effects of the body and 

pressure forces were almost identical when particle material 

density is very much greater than the density of the fluid. 

When both the densities were approximately equal, the effects 

were significantly different. Ragland and Peddieson (19 77) 

extended the work of Healy and Yang (1970) for finite volume 

fraction of the particles and obtained that the loading 

p n 

parameter K (= , p^ and p being the densities of the 

particle and gas phases respectively) had a much greater 
influence on the behaviour of the suspension than the effect 
due to the volume fraction <P . Peddieson (19 76) discussed five 
unsteady dusty gas flow problems in which the motion is induced 
in a semi- inf in it e mass of suspension by excitation at its 
boundaries'. Such flows were shown to exhibit boundary layer 
behaviour. 

In part I of the present chapter, we consider the steady 
gas particulate flow between porous infinite parallel plates. 
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The lower plate is stationary and the upper plate moves with . 
uniform velocity. The lower plate is subjected to constant 
suction (injection) velocity and the upper plate has an equal 
injection (suction) velocity. The analysis indicates that the 
velocity of either phase increases as the suction velocity at 
the lower plate increases. In part II of this chapter, we 
solve the following problem ; 

The upper plate initially moves at a uniform velocity 
while the lower plate is kept stationary. Suddenly, the upper 
plate decelerates or accelerates according to the law varying 
exponentially with time. The lower plate is subjected to 
suction (injection) velocity while there is an equal injection 
(suction) velocity at the upper plate. The problem is to find 
the characteristics of the gas-part ioulate flow in between 
these two plates under such conditions. Closed form solutions 
of the governing equations are obtained by the method of 
separation of variables -and analytical expressions for the 
velocity profiles and skin-friction at the lower plate are 
obtained. When the upper plate accelerates, skin friction 
increases non -linearly with time. On the other hand when the 
upper plate decelerates, skin-friction decreases non~l in early 
with time. For the time t -+ o, the solution for the steady 
case (discussed in part I) comes out as part of the solution. 
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Par I 

Steady gas-pa rt iculat e flow betv ?een por ou s p arall e 1 plates 
3 « 2 Mathematical^ fjor mul .at ion of the P rollem 

We consider the motion of the gas-particulate flow between 
porous Infinite parallel plates, The lower plate at y = 0 is 
at rest and the upper plate at y = h moves with uniform velocity 
u n . The lower plate is subjected to constant suction (injection) 
velocity v w and the upper plate has an equal injection (suction) 
velocity. Equations governing the steady gas-particulate flow 
between infinite parallel plates are the following s 



(3.1) 

(3.2) 

(3.3) 

« 

(3.4) 

(3.5) 

(3.6) 


Here (u,v) and (u , v ) are the components of velocities 

hr Jr 

of the gas and particle phases along and perpendicular to the 
plate length respectively, v and p the viscosity coefficient 
and density for the gas phase, the density for the particle 
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phase and 

18 v 

^ _ -g— 

d p 

sp 

where d is the diameter of a particle and p the material 

sp 

dens it y of t he p art icle s . 


Boundar y_ gp nditj. pns 

(3.7) 

(3.8) 

We introduce the non-dimensional variables as follows 

vt Re 


(i) At y = o, u = 0, " v " = v p =T W ( 00118 "k 0 * 1 '*' ) 

(ii) At y = h, u = u = u m 


- yV'ile’ ~ u 

y = , U = 


u 


, u = -«■ , V 

OO -r CO 


u 


D 


_ d 
- £ 


» p 


p p u ph 

= JL„ p = , Re = -2— - 

P sp p^ ’ ~ y 

P 0 * P P 0 


(3.9) 


where p is the density of the particulate phase at y = h. 
P o 

Using (3.9), eq.s . (3.1) to (3,6) become 


6v _ 

ay 


0 


~ 5u _ dfu 
ay ay 2 


^ p p :(u 


u ) 

p y 


(3.10) 

( 3 . 11 ) 


- 6v _ e v 

V - _ - - 

ay ay 


5 Pp(^ 


V 


(3.12) 


~ ( p n = 0 


ay 


p p 


(3.13) 
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_ QU _ 

v —ja. = f p (u - \L) 
y 6y p 


6v _ 

Y — = P p (y 

p ay 


V 


where 


fT _ 18 

1' — y u. 

Re p 


SP 


Boundary conditions (3.7) and (3.8) become 
(i) At y = o, u = 0, v = v p = v w (constant) 

( ii) At y = h, u = = 1 


(3 .14) 


(3.15) 


(3.16) 

(3.17) 


Here 

h = fRe . 

3 • 3 S olut ion of t h e g o vern ing^ ‘ equ ati on s 

We solve eqs. (3.10) to (3.15) subject to the boundary 
conditions (3.16) and (3.17). Using (3,16) and (3.17), eq.s, 
(3.10), (3.12) and (3.13) give 


? = T p = v w (3.18) 

P p = 1 (3.19) 

Bq. (3.15) is automatically satisfied. 


Using (3,18) and (3.1U), eqs. (3.11) and (3.14) become 



d u __ 9_ u ^ 

ay a;P 




(3.20) 
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du _ 

"v — — -E* =1' p ( u — u ) 

w , ~ v p J 


(3.21) 


6y 


Eliminating from eqs. (3.2o) and (3,21), we get 


3- 


^ 2 - 
6 u 


8 u , /E p ~ \ d~u. ™ Qu _ _ 

“TS r “ O -T9 ~ * (! + p) — - 0 

ay ^ sj ey 


w 


(3.22) 


Eqs. (3.2o) and (3.22) with the boundary conditions (3,16) 
and ( 3 . 1 7) g ive t he f ollowing s olut i on . 


u 


= | (S(l - e Q y/ g cosh |2) + e Q y/» s mh |2 


D 1 + I ) \ + ? ?) 3 coth ir - p h 


P h 


> (3,23) 


u. = 


| {S + e' 


5 Qy/2 GO q]q (y^-Po+Py coth S-^) 

(t N W W 


+ e 


Qy/2 


sinh 


(f p - v^) coth Ph _ 


>(3.24) 
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Coefficient of skin-friction at the lower plate is 

C-f = — (— ) (3*25) 

f YEe dy y=0 


where 


C 


^ g , t 'being the skin— fr ict ion coefficient 


0 pu 
2 00 


or 0 , 


_ 2_ 
RsfRo 


z \ C 1 + =Hf> ooth g 


P 


Ph 


- , P P\ 

v + — t ”) 
w - ' 


(3.26) 


w 


Prom eq. (3.24), we get 


u_(o) = | c 1 + =) € ( 3 - 27 ) 

p 

Eq. (3.27) shows that particulate velocity at the lower 
plate depends on the suction (injection) velocity. 

Por the case of no suction (injection), eqs. (3.10) to 
(3,15) reduce to the following 



u = u p 

(3.28) 

and 

2- 



•"W = 0 

sy 

(3.9B) 

The 

solution of eq. 

(3.29) with boundary conditions (3. 

and (3.17) 

is 



u = i 
h 

(3.30) 


The solution (3,3o) can also he derived as a particular 
case from eq. (3.23) when v w = 0. 
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3 .4 D iscuss ion _of __tjie_ r e suit s 

We present here some characteristic features of the steady 
gas-particulate flow between porous infinite plates, lor 
numerical computations, we use the following values of the 
parameters s 

1 = .8 

v w = ~* 1 ’ ~* 5 

Re =1.0 

lig. 1 gives the variation of the velocities of the gas 
phase and the particle phase with changes in the suction velocity. 
We find that as the suction velocity increases, velocity of 
either phase increases. Also, the velocity of the particle phase 
at the lower plate increases with increase in the suction velocity. 

Part II 

Unste ady flo w__b etwe en p o r ous inf In it e parallel pla tes 
3.5 Matjtema^ical^ fjy m ula tion of _ th e pr oblem 

We consider the unsteady flow between porous infinite 
parallel plates. The upper plate initially moves at a uniform 
velocity while the lower plate is kept stationary. Suddenly, 
the upper plate decelerates or accelerates according to the law 
varying exponentially with time. The lower plate is subjected 
to suction (injection) velocity while there is an equal injection 
(suction) velocity at the upper plate. Equations governing the 
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unsteady flow between porous Infinite plates are 


dv _ 

ay 


0 


~ ~ 2 — 

9u _ h - au _ a^ u 

at 5 y ay 2 



(S - y 


d.v 

at 


+ V 


- av _ a 2 v 


ay ay 




I‘ p (v 


V 


+ 
at 


ay 


( p 
^ r 


v =0 


au - au _ _ _ _ 

+ V - =1 P( U — U ) 

at p ay p 

av _ av — _ 
at p ay L p ; 


(3.31) 

(3.32) 

(3.33) 

(3.34) 

(3.35) 

(3.36) 


Here 


t = 


t u 

o 

h 


and the rest of the variables have same meaning a,s defined in 


(3.9). 

Boundary conditions are the following : 

( i) At y = 0, u = 0, v = v = v (constant) (3.3?) 

jj V'V 

tmm r/Hh 

(ii) At y = h, u = Up = e , P p = 1 (3.38) 

where a is a non-dimensional constant number, a is positive for 
accelerating upper plate and negative for a decelerating upper 
plate . 



66 


3 . 6 S^olutiori p>f_ t he_ oblern 

Eqs . (3.31), (3.33) and. (3.34) subject to the boundary 
conditions (3.37) and (3.38) give 


v = v 


v 


w 


P = 1 

P P 


(3.39) 

(3.40) 


and then eq. (3.36) is automatically satisfied. 

With eqs. (3,39) and (3.4 q), eqs. (3.32) and (5.35) become 


let 


where 


s 4 + T _ f(5 _ u ) 

8t w 8F if P 

(3.41) 

3 u n _ 3u „ 

+ v ■— *- = E p ( u — u ) 

3 1 W 3 y P 

(3.42) 

rR 

0) 

li 

1:3 

(3.43) 

u p = e at H(y) 

(3.44) 

G(h) = H(h) = 1 

(3.45) 

1(0) = 0 

(3.46) 


With eqs. (3.43) and (3.44), eqs. (3.41) and (3.42) become : 
(JH _ V G» - (1 + a) G + I H = 0 (3.47) 

v w H* + (1 p + a) H - F P G - 0 (3.48) 

where * denotes the derivative with respect to y. 

Eliminating H from eqs. (3.47) and (3.48), we get 

(l 3 + a l 2 + b 1 + c) G = 0 (3 .49) 
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where 


let 


D = S 


L, 

ay 


a = (f p + a „ v 2 ) 

v w' 

w' 


~(E + E p + 2a) 


v 


E^ P - (E p + a) (E + a) 


w 


1 = m ~ | (3.50) 

Using (3 .50) , operator with in round bracket in eq., (3.49) is 
transformed to 

3 

m + pm + q 


where 


P 

<L 


ba , 2a 3 
3 2 7“' 


The discriminant A of the equation 

3 

m + pm + q = 0 

is given by 

A = 4p 3 + 27q 2 

or A = 4b 3 + 27c 2 + 4a 3 c ~ a 2 b 2 ~ 18abc 


(3.51) 


(3.52) 


Depending on the values of the parameters, a may be 
p os it ive or negat ive . 

Analyt ically , it is difficult to say about the sign of A 
with varying parameters. Numerical computations have shown 
that positive values of a determine A as negative while negative 
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values of a give A as positive. Tais means that for accelerating 
upper plate, A is negative and for a decelerating upper plate, 

A is p os it ive . 

C ase I 

Let a Toe p os 'it' ive, that is, upper plate is moving with a 
velocity increasing exponentially with time. 

Here, A is negative and then eq.. (3.51) has all the three 
real roots. 



where 
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Hence solution of eg.s. (3.47) and (3.49) using boundary conditions 

(3.45) and (3.46) is 

sinh P;L y x 1 (y-h) 

Gr = ^ e 

sinh y^ h 

a P x y sinh y Pph \(y-h) 

+ b {e ~ e e 


sinh lja • 

sinh v^h-y) i-j_y 
sinh y^ h 


(3.53) 


and 


x i(y-ii) 


H 


s inh y ^h 


n (P-l sinh y^y + Q 1 cosh y^y) 


A ,,P+oc . v w ^1 ^lx 


+ 5 U~rr” + 
a P P 

P i S \y 
e e 


- -i) e 
P 


x 1 h • 

e sinh y^ h 


(P 1 sinh y^y + Q 1 cosh u^y) 


*±7 


sinh y^h 


P^ sinh y 1 (h-y)-Q 1 cosh v^(h-y) 


(3.54) 


where 


— 2 2 
p = h+a + v w X 1 __ X 1 + y l 

1 P P P 


v 


Qp = 


w 


"l 2u i \ 


and 


P P 

A = P(l-P^) sinh y^h - PQ^ cosh y^h 

Pi* 


B = 


3? ( 1-3? ^ ) + p ( ) -~oc g s in h y^ii 


„ p 1 h _ _ 

- Q^P e cosh v^h + Q^P e 


h h 
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Eqs. (3.43) and (3.44), using eqs. (3.53) and (3.54) give u and 
Up respectively. 

At the initial time t = o, the solution given by eqs. (3. 53) 
and (3.54) gives the steady state solution when the upper plate 
moves with the uniform velocity u^. Ibis is exactly the same 
as given by eqs. (3.23) and (3.24) in Part I. , 

Eqs. (3.25) and (3.53) give 



let a be negative so that the upper plate is moving with 
a velocity decreasing exponentially with time. Here a is 
positive and then eq. (3.51) has one real root and two complex 
conjugate roots. 

let pg be the real root and Yg, 6g the complex conjugate 

roots 

P2 = u i + v i 

( U 1 + V 1 i V3 (u^ - yp 

m + h 1 - n> 

5 2 = „ — 2 “ - — g- — - 
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where 


= 


V, 


/C ah a \ /4h +2 7c^+4a'^o-a^b^~18abC',l/2 


168 
2.„_3_ _ 2, 2 


t c ah a u . /4b +2 7c +4a c-a h ~18ahcs 1/ 2 


1/3 
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Corresponding roots of 


D 3 + aD 2 + bD = o = 0 


are 


[3„ = u 1 + y 1 - 


a 


1 ~ S' 

U 1 + V 1 a . J 3 ( u l" v l) 
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y 2 “ ~ g > + 1 { 


/ _ r U ! + T ! a , . 
6 2 ~ { ~ T“^ " 3 } ~ 1 


} — Ar> 1 ]ig 

f3 (u 1 - v ± ) 

{— tA~~ — ■ - > 


X 2 ^ 1 V 2 


The solution of eq.s. (3.47) and (3.49) using boundary 
conditions (3.45) and (3.46) is 

sin v 2 y - x g(h-y) 


& 


sin Vg h 

A ? p ? y sin y ? y p ? h x 2 (y~h) 
+ VT' ( e - ““ — e e 
■“2 sin Vgh 

sin ^g(h-y) x 2 y 


sin ^g h 


(3.56) 


and 


H = 


x 2 (y-^) 


siny h 
2 


(Pg sin y g y + Qg cos h^y) 


+ {~ (F - « + v a - a 2 ) ® 

■°2 P 


P?y 
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P g h x 2 (y-h) 

e 

sin 


(Pg sin u 2 y + Qg cos u 2 y) 


- (Po sijl Vp(h-y) - Q p cos Vo(h-y))> (3.57) 


sin Ugh 

_ - 2 2 
p = + IjlJ , 2 „ X 2 ~ y 2 

2 f f p ' 

„ _ tv p 2 8X 2 "2 


Ag = P (l-Pg) sixi P 2^ “ ^2 ^ C0S y 2^ 


Bg = (P(l~Pg) + p 2 (v w -p 2 ) - a) e d sin v g h 

~ P ? i i _ „ x 2 h 

- Qo P e cos v h + Qo P e 


Eqs * (3.56) and. (3,57), using eqs. (3.43) and. (3.44) give 
respectively the expressions for u and u . We further note that 
as part of this solution, we get steady state solution. 


Eq. (3.25), using eq. (3.56) gives 


°f = 


0 at ~ x 2 h 

2e e 

YSe’ sin Pg h 

2Ag e a * 

B 9 tfRe 


1 ( P o~ X 2^ 

s „ ...... (y e 

' sin Ugh ^ 


+ x p sin Uph -u 9 cos u P h) > 


(3.58) 


3.7 Discussion of the results 


In this section, we present some characteristic features 
of unsteady flow between parallel plates for both the cases. 



73 


For numerical computations, following values of the parameters 
are considered. 

k = .8 



Re = l.o 


Fig, 2 gives the velocity distributions of gas and 
particulate phases at different times when the upper plate moves 
with velocity increasing with time. This graph shows that as 
the time increases, the difference in velocities of the gas and 
particle phases increases. Fig. 3 gives the velocity profiles 
for the gas and particulate phases at different times when the 
upper plate moves with velocity decreasing with time. This 
figure indicates that as the time increases, the difference in 
velocities of both the phases decreases. 

Fig. 6 gives the variation of sk in-f r ict ion at the lower 
plate with time for both the cases. This figure indicates that 
in case the upper plate moves with velocity increasing with 
time, skin-friction increases with time. When the upper plate 
moves with velocity decreasing with time, skin-friction decreases 
with time. 
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FI 6 . 2 VARIATION OF U, Up WITH UPPER PLATE MOVING WITH 
VELOCITY INCREASING WITH TIME. 



VARIATION OF 
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FIG. 4. VARIATION OF Cf WITH TIME. 



CHAPTER 17 


GAS -PARTICULATE PLOW THROUGH TUBES OE 
VARYING CROSS -SECTION 


4 .1 Intr odu ct i on 

Blasius (I 9 I 0 ) analyzed the motion of an incompressible 
fluid through a symmetrical capillary with constriction which 
starts and ends with the same diameter , He developed an 
asymptotic series solution and obtained analytical solutions 
for the zeroth and the first order terms. His analysis 
predicted no extra pressure losses (due to varying cross- • 
section) above the pressure losses predicated by Stokes in 
straight pipes. Tanner (I 966 ) extended the perturbation 
analysis of Blasius (I 9 I 0 ) and obtained closed form expression 
for the pressure loss in a class of viscometers. Lee and 
Rung (1970) investigated the flow in a circular tube of bell 
shaped constriction specified mathematically by a Gaussian 
distribution curve. Using finite difference technique, 
they obtained numerical solutions of equations for stream 
lines, shear 1 stresses, velocity components and pressure losses. 
At Reynolds number equal to fifteen, their numerical procedure 
failed to converge. However, they obtained results up to 
Reynolds number twenty five by using under relaxation method. 
They found the presence of an eddy on the downstream side of 
the constriction when the Reynolds number exceeded 9.9. 
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Man ton (19 71) considered the second order terms neglected "by 
Tanner (I 966 ) in the momentum equations* He developed an 
asymptotic series solution in the ascending powers of S, a 
small positive number which gave the variation of the cross- 
section of the tube with the axial distance. His results 
for shear stress at the wall agreed qualitatively with the 
numerical results of lee and lung (19 70) . lee and lung (19 70) 
and Man ton (19 71) indicated that the point of maximum shear 
stress occured at the upstream side of the constriction. Ohow 
and Soda (1972) considered the flow in a tube with continuous 
constriction or obstruction. They developed the series 
solution in terms of the ascending powers of 6, the ratio of 
the mean radius of the tube to the characteristic length of 
the tube and obtained closed form solutions for first three 
terms of the series. Their results were valid for a case 
when the spread of roughness was large in comparison with the 
mean radius of the tube. In their analysis? they discussed 
seperation and reattachment points for conduits and sinu- 
soidal wall variation. Recently, loorean (1978) and 
Macdonald (19 78) solved the problem in circular tube with 
exponentially varying cross-section of the walls. Macdonald 
(19 78) used Galerkin's Kantrovich scheme and Doorean (19 78) 
developed a series solution by transfo rm i n g the- domain of the 
tube into a rectangular domain. 

Attention had also been given to gas-particulate flows 
through tubes of uniform and varying cross-sections. Kamis, 
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Goldsmith - and Mason (I 966 ) found exp er imen tally that in the 
flow with suspended particles through a tube of uniform 
cross-section, the gas-phase -velocity and the particle phase 
velocity were identical. Later, Karnis, Goldsmith and 
Mason (I 966 ) observed that for large volume fraction of the 
particles, there was a thin layer near the wall where the 
particle density was almost negligible, that is, there 
developed a particle free region near the wall* Sproull 
(I 96 I) observed experimentally that by adding dust particles 
in the turbulent flow through a pipe of uniform cross-section, 
shear stress at the wall decreased, Saffman (1962) carried 
out the stability analysis for dusty gases. His analysis 
showed that the addition of fine dust particles destabilized 
the gas flow, whereas the addition of- coarse dust particles 
stabilized the gas flow. 

Kaimal (19 77) considered the gas-particulate flow through 
a tube with local constriction specified mathematically by 
Gaussian distribution curve. He developed an asymptotic 
series solution in terms of e, giving the variation of the 
cross-section of the tube and calculated first three terms. 

His results were valid for constant gas and particulate 
densities. He obtained identical zeroth order solutions for 
the gas and particle phases. Even at Reynolds number of 
unity he obtained eddy-like motion. Moreover, his analysis 
indicated that by adding coarse, dust particles to the gas flow, 



77 


the shear stress at the wall was increased. This was in 
contradiction to the work of Sproull (I 96 I) . Mathematically, 
he solved an over determinate system of governing equations, 
which probably caused such discripances as stated above. 

In the present analysis, we reconsider the problem solved 
by Kaimal (19 77) and remove some of the short-comings in the 
analysis as stated in the above paragraph. We consider 
variable density of the particle phase and develop an 
asymptotic series solution in terms of small parameter e, 
which gives the variation of the cross-section of the tube. 

We obtain the analytical egressions for the first three 
terms of the series. The present analysis indicates that 
the particles are migrated axially and radially towards the 
wall, thereby resulting in excessive ' accumulation at the wall. 
Our three term series solution does not predict the formation 
of eddy up to Reynolds number equal to eighty. We further 
observe that by adding coarse dust, the shear stress at the 
wall is decreased, which agrees qualitatively with the result 
of Sproull (1961). 

4 » 2 Formulat ion of the pro blem 

Consider the steady gas-particulate flow through an 
axisymmetric tube of varying cross-section. We select the 
cylindrical polar coordinates (R,e,X), where R = 0 is the 
line of symmetry and X = 0 gives the minimum characteristic 
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radius a Q . The radius of the tube varies slowly with the 
axial distance ’X ! viz, 

R = a(X) = a Q S(GX/a o ) , O^G^l (4,1) 

Let (U,V) and (U , Y ) denote the velocity components 

P P 

of the fluid and the particle phases in the axial and the 
radial directions respectively. The governing equations for 
steady motion of a gas~p articulate system are 


UU V + ' 

VUr + .X = 

V (U XX f U RR + B 

: U R'^ 

p 

( u -°p 

) (4.2) 

uv v + 

li. 

R-D 

yy + = 

xl n 

V ( y XX +V RR + R 

V h 

V) 

"p 

(V-V p ) 





R 



(4.3) 

(PU) 2 

+ 

k (OT) 

+ (PV) H = 0 




(4.4) 

Vpx 

+ 

V p V 

= R(u-u p ) 




(4.5) 

U V „ 
P P^ 

4- 

Y pV 

= R(v~v p ) 




(4.6) 

(p u 

p p 

)? 

+ 1 < P p 

V ) + (p Y )_, 
p ' v p p 'R 

= 0 



(4.7) 


where p , p, v are the density, pressure and kinematic 
viscosity respectively for the gas phase; p the density 

ir 

of the particulate cloud and R the interaction parameter 
given by 

tp _ 18 v 

X' — , 

d p 

sp 

Here, y is the viscosity coefficient for the gas phase; 
d the diameter of a particle and p _ the material density of 
the particle phase. ■ In eqs, (4*2) to (4.7), suffixes X and R 
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denote the partial derivatives of the dependent variables 
with respect to X and R respectively. 

I £'P.y. the . fluid phase 

( i) No slip condition at the wall gives 

TJ + a, r y = o (4.8) 

( ii) No flow condition across the -wall gives 

Y - a^-U = 0 (4.9) 

( iii) hue to symmetry of the flow at the axis of the tube, 

Y - 0 and || = 0 at R = 0 (4.1o) 

II Fom the partic le phase. 

( i) No flow across the wall gives 


1 - 0 


(4.11) 


(ii) hue to symmetry of the flow at the axis of the tube, 

au 

Vp = 0 and = 0 at R = 0 (4»12a,b) 

It is to be noted that we have not used the no slip 
condition at the wall for the particle phase. 


4 *3 Method of solut i on 

let us define the stream functions ^ and o for the 

P ’ 

gas and particle phases respectively as follows: 

U | *R • Y = "I +X t 4 * 13 ) 



80 


PU=i* B)P ? 
p p R r pR* M p p 


R *pX 


With eqs. (4.13) and (4.14), eqs . (4.4) and (4.7) 
automatically satisfied. Also, with eqs . (4.13) and 
eqs. (4.2), (4.3), (4.5) and (4.6) he come 


k 

R* 


^RX 

+ A 
R 

S' ^x hs: 


, *RR + T 

r k 

l r 

(♦a 

XX^RRR^ 

- “2 %R + p" 



Rp 

- -£■ 

PR R 

+ 3L 

°R 

^PR 


1 

^X ^XR + 

1 2 

P R 

'XT’ 

? 

h3 *X 

P 

= V 

ri 

Lr 

( ’tec* 

£Rr) " 

r p 


1 

R S 


p 

X PR p X 

1 p 

E V (P P +pax“^,x W + S *pX ( -i ♦pE + p pE^pR 

P p ^pRE / * "p 
1 '1 


- p p W = r P p 3 *R- Jp p 2 "tR 


I ’'pR^pXVVpxR " R *pX ( -~R *pX + PpR* 


n V 


. p * ,„ r ) = ~Rp 3 <|u+ * y 

P p v pRX' p R p pi 


Boundary conditions "become 
( i) At R = a(X) , 


1>=%, ^=0 and * p ***p 


o 


(4.14) 

are 

(4.14), 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


( 4 .19 a, h , c ) 
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( ii) At R = 0, 

R = °» (r *r)r = 0 ( ^ *20a,To, c) 

% = 0? R X = ° 5 (r p“‘ ^p3.^R = 0 -(4*21a,b,c) 

.P 

From eq. (4,1), we notice that when e - 0, the tube 
is of constant radius a Q provided S(o) = 1. lor a non- 
zero G , the variation of S in the axial direction depends 
on GX instead of X alone. 

We nondimensionalize the variables as follows: 


r _ GX R 

r M He jm-4 

a o a 0 

~p ” P p/ P p 0 ’ p 


H 5» n> = 'l>/'P 0 3 $p($, n) 


^pAp o ’ 




i a .p 


(4.22) 



being the constant density of the particle phase. 


Using (4.22) in eqs. (4.15) to (4.18), we get the 
governing equations of motion in the non dimensional form 
as follows: 


Vnr rz + % > + 6l li 


_ e 
2 


, 1 , 1 . ^ 1 . a - . , 

d) + — d> — "-ft- 6 + 9 — — p T •“ T nn 

n n^wm $ m n 3 n n p n n -P n 


I +d„ (4-23) 


E e^Z***M ~ VcU 3 +{> 


Ur 


P P 

4> ririr + - 4> 
n ts.9 n grin n 


6 2 ~ 6 In ^p^S* e ri ^Pl 


(4.24) 
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~n ^Pn ^Pn? P p rT $ Pn P p? + “ ^P£ %>n p Pn + JJ ^Pn P p 


e - - -3 

n ^pnn P p P p v n 


Pn ^ ~ P p P ^ 


(4.25) 


,2 - 


e 2 2 


~ * Pri % P pf V ‘k % “ T *pn *p« P p " ;~2 * P S P p 


_2 


_3 


"n P P *PCn'*P5 6 ^ P P ^PC ~ 6 a P P \ 


(4.26) 


Here 


H = — , .Reynolds' number of the flow 

cl o 


P 


a = 


— 9 

0 


... P 


Pb 

- ZlSL 

V 

9 a 

~ " S 
P 

P 

% p 

Pa 3 

0 

* 'o 

2 

r\ 

Pa 3 

0 



P Pn’ 

P 

" *rT f 


Boundary conditions (4*19) to (4,21) become 


At n = s, 

4> = lj ^ = 0 , 4>p = ^ ' (4.27) 

At n = 0 } 

<f> = 0, ~ * 5 = 0, (~<f> n ) n = 0 (4.28a,b, c) 

* p = o, i * pc = 0, (i i tp ) n = 0 (4.29a,b,o) 

P P 

4 . 3 Series expan sion 

In the present analysis, we take R e to be small so that 
SR = 0(1) . This gives low Reynolds number solution in which 
viscous forces dominate the non-linear inertial forces. 
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Here, we develop series solution in terms of the 
ascending powers of e and calculate equations governing 
the first three terms. Since e -* 0 gives the tube of 


constant radius, we expand functions 4 , <j> , p and q 

P P 

asymptotic power series in e as follows; 

in 

<*) = + e 

o' 1 ) +e 2 ,(2) + .... 

(4.30) 

. _ .(0) , p ,(1) , p 2 , (2) 

♦p “ *p G(i> p + 6 . *p + 

(4.31) 

p = p"(°) + g 

p p 

+ E 2 -(2) + .... 

P P 

(4.32) 

1=1 3?(0 

+ P j n ,G ) 

(4.33) 

where 

q(°) + Gq.^ + e 2 q^ 2 ^ + **•• 

(4.34) 


Since the zeroth order term gives the classical 
Poiseuille flow and for such a flow, the pressure must he 
a linear function of the axial coordinate. Also, in 
dimensional variable first term in (4.33) becomes XE'(GS) 
which has a constant gradient with respect to the variable X. 
Hence non-dimensional pressure q can be expanded in the 
form of (4 ,33) . 

Substituting eq_s. (4 .30) to (4.34) in eqs. (4.23) to (4.26) 
and collecting like powers of 6, we get the following system 
of differential equations governing various order terms. 
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g f i o) ^ 4 o) p‘ o) -^ i o) } 

1 _(o) (o ) , (0) ,(0) \ 1 (0) f .(0) p(0)„ (o) p(0)) 

~~2 Pp 4 ^4 " ^de - n*P5 C *Pn PS *PS % ' 


'PS 


s p5n 


’PS 


+ .i p( 0 ) *( 0 ) *( 0 ) 
p P P W 9 Pri 


S e con d ord er equ ations 


1.2 *i !l 


(4.48) 


- 2 *(o) -(8) = 4 (i) + s o^ 1 ) i 1 ') 

t) v n P p p n P 


, R r 1 . ( x( 0 ) a^ 1 ) + 

+ fi e L “2 1 ^ 9 n 

♦<°> ) - % (♦<<» 

«(!) + 
v ri ri 

4 1 ) c (0) 

S pn 

+ ±- 
p , - ) 

(/°> + 4 

, ( { 1} 4 0) >> 

. i 

P 

4 ( 0 ) 

v «» 

(4.43) 

a ( 2 ) = _ i_ 
4 " „ 2 

(" 4 in - 4V 

~noc p£ 0) ^ (1) ) 



-i( 

n 

p 4 0) 

o (1 )) 

P 




+ ?e 

p 3 

(p 4° 5 4°V - 

,(°) a(0) 

n * 5t > % 


(4.44) 

«- ^ o)S e } 

- F ?<o)%< 2 ) + 

H P p Pp 

p(°)(3a 
p v P 

4,(0). 

-BF^) 

-( 2 ) 

P 

= - 3a 

P 

p p 0) ' 

M P> 





+ F -4 X) ( 

4 ( 0 ) p^ 1 )+ 2 « 
pp P 

(i) 

Pp 

F (0) ) 

P 



i m (o) ~(o) A (o) a^ 1 ) p^°V 4°) 4>^°r p^ 1 ^) 
+ ^ ( 4n W P p + > %* P P fc» W P 
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p(°)+$(°) p(°)- 

■Tk V r y~\ T ^ r r\ T n ^ Tn_ 'V^ 


pn Pn 


Pn 


Pn 


,(0) (0) 7 (1) } 

T>n Pn P5 ; 


+ 1 uWJo) -(o) + ,(o) (1) -(0) + 4 (0) (0) -(1)) 

+ n V *pn p pn + *p 5 *Pn p Pn + *PC ♦ n Pn " ' 


pn 


* v*™#’ r* *s, 11 c f- 4 ? 4 °>> 


( 4 #45 ) 


a p(°>V>- fT p(°) 2 / 8 ) + p(°)(3a B 2 > 

P ? P PC P 5 P^ P 


■3a p(0) ^ 0) + 4 0) f ))+ ^ 1)( i 0) 'SM^p 01 ' 


+ 1 (,(1)*(0) j(0) + ,(0) (1) -(0) (0) (0) 

n U Pn PCS P P" P5? P Pn P 55 

i c <5 * Co) p(0) + (j, ( o) ( i) ~(o) + (o) (o) 

" n { Pn PC pn $ PC P PS *Pn ? PS 




p- ( P) 

H pC 1 


i ,,(i),(o) -Co), A (0).(1) -Co) + ,(o) (o) 

(Clnr P nr/ + *.DF *DE P pn *P 5 *P C 


PC 


p'~'+ A' '■''A 

! pn 9 PC S PC 


p( x >) 

*pn ' 


^ 1 ,.(1) (0) r(0), 6 (0),(1) -(o) .(0) (0) 

+ 7 Op 5 +p?n P p *p5 0 p?n’ p P *P? *P5n 


;(!) 

P 


pO ') 


1 / ( i) (o) -(o), .Co) .( 1 ) ~(o) + Co) (o) 
^2^ p 5 ^p g ^ *P C *P C P P ^ 5 


P (1) ) 


Boundary conditions (4.27) to (4.29) become 


(4.46) 


At 


- S. 


k( n ) = o for n = 0,1,2, ••• 
n 

^°) = l, <}>( n ) = 0 for n = 1,2, ... 

(°) = i 4 n ) = 0 for n = 1,2,... 

'p P 


(4 .47a, b, c) 



b Y 


At n = 0, 

,0) 


kO) = 
p 


°> $ Y } 


o, - 1 - 4 (n) 

n v p? 


o, (i e>),=o 

for d = 0 , 1 , 2 , .. . 
0 f or li = o ,1, 2, . , . 


(4.48) 


.Boundary condition for various order terms in symmetry 
condition (4.29c) at the axis of the tube “become 


(•i _J- T a(°)) = o 

^ ~lcf) V Y u 

P 


(4.49) 


["A 1 . „i_ (p(0) 1 ) _ *p( 1 ) 4> C 0 ) ) -] _ Q (4.50) 

L ’ ^ -AM 2 P Pn P P -TV ' “Y V ' 

P P 

(i_i r?(°>V 8 > - p(°)(-p( 2 )(,(o) + pC 1 )^ 1 )) 

n /^n3 u d Dn P v P On P Pn ' 


n p(o) 

P 


P Pn P ' P Pn P Pn 

+ p^ 1 ) "3> = o 

K n n tv --■ * y~i 


(4,51) 


4 .4 Sol utio n for the ^ var lo r de r^ term s 

( i) S olution of zeroth, orde r equations 

Using eq_. (4.37) in (4.35) and then differentiating 
with respect to n , we get 


" 3 ^nnnn ~ ^ + 3,1 “ 3 *( 0) = ° 


(4.52) 


Solution of eq.. (4.52), satisfying the boundary condition 
( 4 »47a,h) and ( 4 .48 a, b , c ) is 
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4 > 


( 0 ) 


“T 


s 


s 


(4,53) 


In case, we solve eqs. (4.37) and (4.38), we find that 
zeroth order density becomes unbounded at the wall. Hence, 


we assume that 


= constant. Karnis, Goldsmith and 


Mason (I 966 ) found experimentally that in a tube of constant 
cross-section, the density distribution remained uniform. 
Our zeroth order solution represents the flow in a tube of 


constant cross-section. So, the assumption pI (JJ = constant 


~(o) 

P 


agrees with the findings of Kamis, Goldsmith and Mason (1966) . 
In this case, eqs. (4.37) and (4.38) become linearly dependent 


and the solution of is 

(0)_ n 4 

♦p “ ”2" ~ 


(4.54) 


and 

= j5*/a = const. (4,55) 

Eqs. (4.53) and (4,54) show that the zeroth order stream 
functions for gas and particulate phases are identical and 
they represent the Poiseuille flow when S is a constant. 

This is in agreement with the experimental results -obtained by 
Kamis, Goldsmith and Mason (19 66). 

Eq. (4.36) shows that q ^°- > is a function of % only. 

With eqs. (4.33) to (4.55), eq. (4.35) gives 

( £ E(C)) 5 = - (4.56) 

s 
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( i i ) SoiU rtj.q n of the f_^jt_order_ equ ations 

With e qs . (4.41) and (4.53) to (4.55), eq. (4.39) 

he comes 


--•“’-■m ; 1 * js • ( , i> ■ .?> 


<t> 

n rinn 


o(0) o 4 

+ 38B B ( 1+ .J2 r -)S»(Sf-^-- 3 =g) 

P s 1 S J s° 

(4.57) 

Differentiating eq, (4.57) with respect to n and 
using eq. (4.36), we get 

-3 _ 2 n 2 *( 1 ) +n /l) _»( 1 ) 

11 " M nnn Y nn n 

?t (0 ) 5 7 

= 32R ( 1+ -2 — ) S' (3U - ^-) (4,58) 

P . S S* 

Solution of eq. (4.58) with the boundary conditions 
(4.47a,b) and (4.48a,b,c) is 

/ 1 \ qi 4 r )2 n^= o r)6 -i n® 

* = e 6 f-< 1+ 4 r~) <9 p + ? pr - 9 jir) t 4 - 59 ) 

Eliminating 5U) from e qs . (4.41) and (4.42) and 
substituting eqs. (4,53) to (4.55) in the inhomogeneous term 
of the resulting equation, we get 

' *(u + 41 ,d) = 16 an « 4 0+ 4 + 4) 

P £ S DH r< 3 .0 D oh o D cl 4 




S' 


/. 5 p°\ /S' a S"w4 n 2 « 4 2 " 6 It 8 , 

- R e (l+ "t- )( “ s ^ " ”s” 9 + > 

( 4 *60 ) 
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Let S be a linear function of 5 so that S M = 0 
Eq . (4.6o) simplifies to 


AA , n Sl A 1 ")- v ( -i , p i°\ S' 2 ,4 n 2 n 4 ,2 n 6 

p 


PC 


s 9 g2 ” g 4 + 3 ^ 


q n8 
9 


) 


(4.61) 


Using Monge’s method, solution of the eq. (4.61) is 

x(l)_ R s' (1 . j°\/4 T ' 2 u 4 , 2 n 6 1 A + 

*p - R e <T (1+ ^ )( 9 + f( s } 

(4.62) 

wh-ere f is an arbitrary function to he determined from 
the boundary conditions 

f(0) = f(l) = 0, f(o) =0 (4.63) 

With (4.47c), (4 .481, e) and (4.63), eq. (4.62) gives 


(1) S' p p°\/4 ^ A 2 A in 8 
*p R e g~ < 1+ ~ ~5 + 3 ^ “ 9 

r|2 

+ ~yr *-> —tf 

s 

Using eqs. (4.53) to (4.55), (4.59) and (4.64), 
eq. (4.41) gives 

,2 K * r ,„2 _2 


(4.64) 


p(D= _ I®- (1 _ » j _ L (i- Sg_)ca- -V 

P p-S 3 S 2 2a S 2 S 2 


-1 


(4.65) 


The p^ 1 ), in eq. (4.65) is infinite at n = S, Indicating 
P 

that the particles accumulate on the surface. 
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Using eqs. (4.53) to (4.55), eq. (4.4o) gives 


,(!) 


-32 n 


S' 


(4.66) 


,(0) 


q'g-' is known from eqs. (4.57) and (4.59). 

( i i i ) So l ut ion f_or _ seoon jd or de r_ equ ations 

With eqs. (4.45), (4.53) to (4.55), (4.59), (4.64) and 
(4.65), eq. (4.43) becomes 


i *L 8 > - % ^ *< 2) = 4 ( U~ ss&wo 


n ' nnn 


,3 


£. 

, 2 - 


+16 Q R e -401^(1^) 3 D 

D b b b o 

7r( 0 ) O O A 6 

.„2 Mi V S'*, 40 212 ^ , 114 2q8 

_ 4Re(1+ .i_) + 

, 38 r) 8 

+ r pi 


(4.67) 


where 


iP, 


o F p 




Differentiating eq. (4.67) with respect to n and using 
eq. (4.66), we get 


l3 - s* 2 ^ (2) + - S 4 2) 

= 325 H e §^LU0^U-^) 2 -(3-^'J 
+Re 2 (l+ -l—) 2 4 -152 ~g + ^4 - "g) 


5 1 2 

+ 320 ~§- ^5 


(4.68) 
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Solution of eq. (4.68) subject to the boundary 
conditions (4.47a,b) and (4 .48a,b, c) is 


* (2) =§ 


2 n , 

S' ~s (1 - -j) 


P 2 2 


+ Q R 


+ 20 


S 

e S’ 
n 4 

,T‘ 


S' 

, p(°) 2 


24 


n6 


so n " . ” 8 - 1 v 

r + Y 


2 1,2 1 n 4 In 6 -, 

3 ^ + ? 7 * ~ 9 ^ J 


p( 0 ^ p _ t 2 y<2 oio 

■R^ (1+ ) §U.— ’ ( + 

P'S 2 S 2 12700 540 s 2 54 s 4 


+R‘ 


479 n2 . 53 n4 


19 




13 


n8 


19 


nlQ 


" 3? ' S T + 90 - IS® po (4.69) 

With eqs. (4.41), (4.53) to (4.55), (4.59), (4.64) 
and (4.65), eq. (4.45) becomes 

*(0) p(0) 2 ^ 2) = -ns ^ 0)2 (4 2) t£ 0) + ^ ^ X) )-F 


a 


3 , 2 ± 


p 2 3 


+16 Si SE (i- 2-f + ?£) - 640 -V i (1- 

g5 ~ g2 g^= g^ pf g^ 


— 4 R . ( 1 


p(°) 


.)(J5. _ 212 f ', + 114 


T)4: n6 


a 


9 g 2 ' "3 s 4 

n8 


208 n 
"9 ^ 


38 S * 2 

jr 


(4,70) 

With eqs . (4.42), (4.53) to (4.55), (4.59), (4.64) and 
(4.65), eq. (4.46) becomes 



a 


.(°) ^ f L ~ C 5 ?< 0)8 c*< 8 > tfhj?) 4 0)S i 2) j 


S' 3 n 2 
|3S ' S 


+ 188 It? t? (1 “ p? )(1+ fsr - ;i : > 


2 8n 4 


S 


6 S * 2 ^ 2 rft 

~ 8 " ~T* :2 ( 14 -58 -g + 52 ~|) 

cd S 


a S' 


4 | (1+ !eZ )( 16 ^ _ 9|8 n* + 588 _ 138 

a S° e o s 2 9 r 4 3 s® 


o(0) 


+ 


,8 

"8 


, 256 ti 10 v 
+ ™ ) 


S^ 

.(4 .71) 


r( 2 ) 


Eliminating from eqs . (4.70) and (4.71), we get 


6( 2 ) + T§l . ( 2) 

V S *p n 


16 


9 2 

S' 1 ^ 


n 2 n 4 


hr* 7 p?B 6 ~ 26 p? + 24 ^ 4 ) 


S 

1 S' 


+ 8 i s “3 7- J)(4-15^ + 13 2^)3 


. „ n S' 2 " S r- ,5 2 n 2 1 n 4 1 i 6 
+ E e Q ^ 1-^18 " 3 gT + 2 “? ‘ 9 ^ 


.41 §1 (|9 _ g4 ; 
S 3 


l2 . 2n " 4 80 " 6 . f\ -l 

+ 20 ^ - t ~e 1S> -J 


+ 


p ( 0 ) qf 2 _2 - 2 ri 2 -» -i n 2 

E e( 1+ “^)Di8 - 9 ^r) 


+ 


8 S 1 I M n 2 2 n 

9 g 5 Y 


s t s 

4 SI I ,104 _ 710 fi 474 jf _ 1034 1. + 218 

" 4 g 3 B 9 9 p2 + 3 g4 9 g^ 9 ^3'* 

9. . ^p°\2 S' 3 n2 1 818 479 ^ ^ 53 n4 

- R e (!+ -Jj—) jr ps < 2 To 0 " F 4 o + 54 
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Using Mange* s method, solution of eq. (4.72) is 


h s >= ♦<*) _ | s - s ( i-4' 2 

p 0 s 


s<°> 


S' n 

"e v ~ ‘ - > g~ 3 

P b 


) + R (1+ -2— ) 


xf T SJ /l04 710 n2 , 474 n 4 l 0 34 i 6 218 n 8 

* S 3 ~9 Tp + ~3~ ^-~9—^ + T s 


_ ( i- ^ (4 - ig) n4 _ i ^ + § i * (1 _ 2 > 

- £ ? fi J CW J ^ ^ f a- J 

2 4 

(4-15 ~?r + 13 ~2|) H+f^( T1 /S) (4,73) 

where f ^ is an arhitrary function to he determined from the 
boundary conditions. 


f l(0) - f 1 (0) - 0, %(!) ”3 "^3 

P s 

* # • 

f l f l 

— *■ — ► 0 as n -* 0 

ns „2 


conditions (4,74) give 


n 4 n 2 


f l (n/ S ) ^ ^ '? + ? 


(4.74) 


(4.75) 


4>( 2 ) is now completely known from the eqs. (4.73) 
P 

and (4.75) , 

Using eqs • (4.53), (4,59) and (4.69), eq.* (4,7o) 


gives 
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;( 2 ) 


1 r ^ S ' 2 „ ,, P i°\, 

' { 4 :T“ k e ( 1+ —“) ( 


( 1- 


96 2240 7560 ^ 

W" ~ 9*' c 2” + 9~“ 0 4 


7680 n 6 -2112 ti 8 v 
-9™ ^ + -F- “8') 


+ 36 fz R e U+ ”t~) nz < 4 - fz) 
S P S s 


p(0) 


- . 2 n 2 n 4 

3, y ^ ( 1-4 -£ + 3 ) 

b b 


5 « ' 

' 6 ' 


-41 ( 3 - |- 5 - + 46 \) 


16 S' 2 ,38 

Y g ^ 3 
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„ i§i n , 258 n 4 _ 134 n b v > 


^ TE " "T" 
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rfi n 2 
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(4.76) 


Using eq. (4.69), eq. (4.67) gives 


4 D = - g (|8 -80 4 


S' 


*Pel 

s s 


2 p ( °) 

, 1240 S'* v 2 M V . 

+ ^ R e (1+ -^~) 

Substituting eqs. (4.53) to (4.56), (4.59), (4.64), 

(4.65), (4.69), (4.73) and (4.75) in eqs. (4.30) to (4.32), 

_ 2 

we get solutions of $,■<(> and correct upto S . 
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Shear stress at the wall 


Shear stress at the wall 1 = a (X) is 

V + ^ a jffi" a XX^ )a s3/( 1+a |^ 


(4.77) 


where ^R 821 4 are the stress components given by 

a nn = « P + 2PV PJ 


P + 2pv 


a = pv ( 3g + 6J) 

xr ^eR ax ; 


Conditions in (4.27a,b) and (4.28a, b,c) give 
At R = a(X) 


V = 0 » W ®X» h 


“ a £ %R 


(4.78) 


Using (4.78), eq. (4.77) becomes 


V DV( S ^XX ~ ^ % + I ^RR^ 


(4.79) 


Ron -dimensional stress r at the wall is given by 
3 m 

x = i°.i. 

W lJ) Q PV 


% = k ~ ^ + I- ♦« 

p(O) 

r w = - % + |~ e R e (! + -^-) K 


-e 2 C%(|s« 2 -l) + 4^QR e (l2T^. + ^) 

S O u . 


Q *(0 ) P 

+ tt! R?( 1+ ~?r~“) 2 ^h“ U+o(e J ) (4.8o) 
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Ve locitx P rof_ile _s 


The axial velocity U, for the fluid phase is given hy 


- i 

U = - x - 
n 


J n 


3 (1 "* ^ + e R e % ( 1+ H“) 
b .b S 0 


r> 4 3 


+ s 2 l:^ s ' 2 


Q n ^ 

X (*>« 4: — w* 4- 4: - — — 

^ S 4 9 S 

i2 n4 




r’ 1 ~ 4 rs + + Re q - 3 

O D S S 


q 1 Ko n2 n4 n6 f|8 

x (t ( 5 ° _ 96 ^ + 120 -160-^ + 10 ~g) 

s s s s s 


. Y / lO 8 . rf* 3 

- Y( 18 “ 3 ^ + 3 3 " Q 

n(0) 


8 n 

s 


2 


A r* fjS 




rTv ^ 7 9 9 A 

J p 2 Mj _ P p / 1636 19 16 ^ , 318 n 

+ H e ( 1 + -K-~) -y [yfffi - — £5 + S 5 ~ T[ 


152 ^ 


130 


3 T 90~ ^ " B‘50‘ s 


" 8 228 ' 1 ^)3+o(e s ) (4.81) 


The axial velocity U , for the particle phase is 


given hy 


- _ 1 1 

P “ tT 


"pn 

,2 2 T” 2 In Q! 2 


V i. + ^(2 2* -i) + e 2 t ~ f ? 

P p 


, . / 64 IS' _ 64 T 21 H + 

+ 4r(g ,W5 ‘ y* ‘"“'Tj 2 Y gO 


S" w s~ S 

2 64 IS* 


n 2 n 4 
+30 ~v) 
S 4 
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((3-86 4 +36 4) + f a; (8-76 4 + 168 ^ ■*>* %) 

8 S . S S S S 


1 ci n2 n 4 

18 & e (l+ ^-~-)(8~52 — g + 66 - 

S p S 

2 - -(0) 


7^6 

-16 -~g) 

s i S D 


^ n'- u V O 6 

-1- SJ 2 p / -i , n p wl92 2736 ^ 8352 nD 

h s r ~y iL e 1+ *- 4 r-)(- 9 »" - - 9— ^ + -9- ~s 


§160. ^ + .§p, ^ -) > + 0^3) (4.82) 


4 • 5 S.isjili. 3 ision^ of ^ t he re sults 


In this section, we present some basic features of 
the gas-particulate flow through a tube of varying cross- 
section. lor the illustration of the results, the following 
form of the tube is considered. 


S = 1 + .2? 


lor numerical computation, the three forms of the 
series expansion (4.3o) to (4.32) are considered with the 
following values of the parameters 

R = 1- 125 
0 

1 = 1.0 
Y = 1.0 
6 = 0 .1 
Q = .25 

fig. 1 shows the stream line pattern for gas and 
particle phases at R @ W 1.0. She figure indicates that the 
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stream lines are almost parallel to the surface. This 
figure further indicates that stream lines for the particle 
phase lie above the stream lines for the gas phase and as 
we move towards the wall, stream lines for both the phases 
approach closer and closer to each other. Moreover, 
the flow exhibits no eddy like motion at this Reynolds 
number* Rig. 2 gives curves of constant density in a 
meridian plane. This figure indicates that dust particles 
have tendency to move towards the wall, which results in 
excessive accumulation of the particles there. Due to the 
curvature of the pipe, centrifugal forces acting on the 
particles causes the particles to move towards the wall. 

In fact, due to the accumulation of the particles on the 
wall, particulate density becomes infinite there and the 
present analysis fails* Rigs* 3 and 4 give the stream lines 
pattern for the two phases and constant density curves at 
R = lo.O* The flow variables follow the same trend as 
for the case R = 1.0. Rig. 3 indicates that there are 
no eddies in the flow even at this value of Reynolds number. 

Rig. 5 gives the axial velocity distribution at 
C= 1.0, 2.5 and 5.0. The figure indicates that > U 
at all the sections of the tube, that is, particles move faster 
than the fluid. This figure further indicates that the 
difference in velocities for the gas and the particle phases 
decreases as we move downstream. As expected, velocities 
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of "both the phases decrease in the downstream direction. 

It should he noted that for the particle phase, we have not 
improved the 'no slip' condition at the wall of the tube. 
Computation of the result show that the particulate 
velocity at the wall is almost equal to zero. 

Pig. 6 gives the shear stress distribution for the 
gas phase and particle phase. The figure shows that shear 
stress decreases sharply in the initial portion of the tube 
and then, it decreases slowly in the downstream direction. 

This figure further shows that the shear stress for the 
gas phase is higher than the shear stress in gas-particulate 
system, that is, addition of dust particles de creases the 
shear stress at the wall. This result is in agreement with 
the results of Sproull (1961) and Saffman (1962), However, 
the predicted decrease in shear stress is less' as compared 
to the observed decrease in shear stress by Sproull (I 96 I) 
for the same amount of dust. This may be due to the reason 
that Sproull considered the turbulent flow, while the present 
analysis is valid for laminar flo?? only. In the present 
analysis also, Pig. 6 indicates that as the Reynolds number 
increases, the percentage of decrease in shear stress increases, 
further, decrease in shear stress by the addition of particles 
implies that a given flow rate can be achieved for a gas 
with dust with a lesser pressure gradient than the pressure 
gradient needed for clean gas. This result is essentially 



similar to the corresponding result of gas-particulate 
flow in a pipe of uniform cross-section. 

Since the leading term in. the present analysis 
represents the "balance "between the viscous and pressure 
forces, its solution, as expected, is the Poiseuille flow. 
Effect of the inertia terms appears as perturbation in the 
higher order terms, Thus the analysis from this 
consideration is strictly valid for low Reynolds number. 

In the development of the analysis of the present investigation 
there is no restriction on the magnitude of Reynolds number. 

Eor higher Reynolds numbers, the contribution due to second 
order terms is larger than the contribution due to zeroth 
order and first order terms which results in formation of 
eddies near the upper wall of the inlet portion at R e =125 . 

Eig. 7 clearly indicates the formation of eddies. In 
Eig. 8, we have drawn the velocity profiles for gas and 
particle phases. Eig. 8a and 8b indicate the reverse flow 
near the wall for gas and particle phases respectively. 
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C HAPTE R 7 

CAS -PARTICULATE PLOW THROUGH A CURVED PIPE 


5 .1 IntiAoUu^ijDii 

In a cross-section of a curved pipe, under the action 

of centrifugal forces and pressure gradient, a secondary 

motion is set up, which greatly reduces the volume flow 

rate • Hence to get a certain volume flow rate, a larger 

pressure gradient is required in a curved pipe than it is 

needed in the straight one, Eustice (1908) first 

investigated experimentally the presence of secondary 

motion in the flow of a curved pipe. Dean (1927) solved 

the E'avier ~ Stokes equations governing the flow in a curved 

pipe and ohtained the solution for the first approximation 

2 

only ( i «e . terms of order A? are neglected where a is the 

R" 

radius of the' cross-section of the pipe and R is the radius 
of curvature of the pipe). His analysis also indicated the 
presence of secondary motion and the results were in 
qualitative agreement with the experimental work of Eustice 
(1908). Later, Dean (1928) considered higher order 
approximations and found that for a given pressure gradient, 
the volume flow rate depends upon the curvature of the pipe. 
Here, Dean (1928) did not consider the effect of secondary 
flow on the volume flow rate . 
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Dean (1959) derived analytical expression for the 
volume flow rate by considering the secondary flow through 
the curved .pipe of a circular and a rectangular cross- 
section. Here, Dean (1959) replaced the secondary flow 
by a uniform stream and found that for a given pressure 
gradient, the secondary motion decreased the volume flow 
rate. He further observed that the point of maximum 
velocity was shifted outward. Smith (1975) obtained an 
asymptotic series solution of Navier -Stokes equations for 
large values of Dean number and for a fully developed 
laminar flow through the curved pipe of a circular, a 
rectangular and a triangular cross-section. All the above 
authors assumed that the radius of curvature of the pipe 
is large in comparison to the radius of the cross-section. 

Hawthorne (1951) analyzed the flow through a bent 
circular pipe. He showed that secondary circulation 
remained unchanged if stream lines were geodesic surfaces 
of constant total pressure. He further observed that the 
secondary flow in bends was oscillatory. 

To the author’s knowledge, no attempt has so far been 
made to study the nature of the two phase flow through a 
curved pipe. In the present analysis, we extend the work 
of Dean (1959) to gas-particulate flow' through a curved 
pipe for studying the effect of secondary motion an the 
volume flow rate at a given pressure gradient. We assume 



that (i) the radius of the circle in which the central line 
of the pipe is coiled is large in comparison with the radius 
of the cross-section of the pipe (ii) the actual secondary 
motion of the gas phase is a uniform stream in the central- 
region. we calculate the stream function for the particle 
phase from the governing equations. Due to the presence of 
the interaction term in the momentum equation for the gas 
phase, the method of separation of variables as adopted by 
Dean (1959) failed to give closed form solution. As such, 
we solve the governing equation by Green's function approach 
and obtain closed form solutions for the pipe of a circular 
and a rectangular cross-sections. Dean's (1959) case comes 
out as a particular case for the clean gas. 

We find that for a given pressure gradient, the secondary 

flow effects the volume flow rate in the presence of the 

particles also. When the particles are fine in nature, the 

volume flow rate for the clean gas is more than the volume 

♦ 

flow rate for the gas-particulate system. For coarse 
particles, the volume flow rate for the gas-particulate 
system is more than the volume flow rate for the gas flow. 

In both the cases, it is obtained that the region where the 
primary motion is greatest, shifts outward. 

5.2 Mathema tical f o rmulation of the problem 

Consider the steady gas-particulate flow through an 

She radius of curvature R 


infinitely long curved pipe 


of the pipe at any point is large in comparison to the radius 
of the cross-section a. fig* la shows the system of coordinates 
that has been found convenient in considering the flow through 
a curved pipe. OZ is the axis of the anchor ring formed by 
the pipe wall, G is the center of the cross-section of the 
pipe by a plane through 02, that makes an angle $ with a 
fixed plane and GO is perpendicular to OZ and is of length B. . 

If P is any point in the cross-section and the distance 
of the point P from OZ is r, then the position of the 
point P is specified by the cylindrical coordinates 
(r,8,z) * The surface of the pipe is given by r = a, a 
being the radius of any section. 


Here, all the flow variables except the pressure are 
assumed to be independent of 4 1 . If (U,V,F) and (TJ ,V ,W ) 

Jr hr ir 

are the components of velocities of the. gas and particle 
phases respectively, then the governing equations of motion 
in cylindrical polar coordinates for a steady gas-particulate 
flow are 
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(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 


Sr p p 

Here, y , p and P represent the kinematic viscosity, 
density and pressure for the gas phase respectively; p p , 
P g p and d the particulate density, material density and 
diameter of the spherical particles respectively and P 
the interaction parameter given hy 


P 


18 y 
“TT" 
d p 


sp 


y being the viscosity coefficient for the gas phase. 

The radius of curvature R of the pipe is large as 
compared to the radius of the cross-section. As a consequence, 
it can be considered that 



V/Ith these approximations, eqs. (5*1) to (5.8) 
to the following form: 
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Prom eqs. (5.9) to (5.11), we get 


reduce 

(5.9) 

) (5.10) 

(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 


P =-AR<J> + g(r, z) 


(5.17) 
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where 

A = ~ R 5 T ^ 0 (5.18) 

A gives the space rate of decrease in pressure along 
the central line of the pipe and comes out to he a constant 
because R > a. Rue to this assumption, the factor 
~ in the term of equation (5.1o) is replaced by g* • 

It is stated by Eustice (1908) that to cause a given 
rate of flow, a larger pressure gradient is required in a 
curved pipe than in a straight one, the difference being 
considerable even when the curvature is small. Here, there 
is not a constant pressure gradient as there is in the case 
of flow through a straight pipe. But it is natural to 
define the mean pressure gradient as the spaoe rate of 
decrease in pressure along the central line of the pipe. 


In a curved pipe,, the primary motion is along the line 
of the pipe and the secondary motion is in the plane of the 
cross-section. In the present analysis, we calculate the 
flow rate for the secondary motion. 


We define below the stream functions 'f* and ip ' , for 
the secondary flow for the gas and particle phases respectively. 

U = - * W = ( 5 .19 a, b) 


a z 


5r 


ai> 


3 ^ v 

p U = - ~~ E. . P W = 

p p 02 ’ p p 3r 


( 5.20a, b) 
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•*"ith eqs . (5.19) and (5.20), e qs . (5.12) and (5.16) are 
automatically satisfied. 


i + (2L.4- . 2* S_)v + — E,(T_T ) (5.81) 
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Eliminating P from eqs. (5.9) and (5.11) by cross- 
differentiation and then using eqs. (5.19) and (5*2o), 
we get 
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with e q s. (5.19) and (5.2o), eqs. (5.13) to (5.15) become 
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(5.24) 


pOr ~ p p dr 


(5.25) 


In the approximate equations (5.21) to (5.25), the 
coordinates r and z appear only in the operators 
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8 

6*z 5 hence any point in the plane that is selected 
for ref ei once can he taken as the origin, and jib denote 
X,T for the corresponding coordinates (fig. lb) . 

he n on-dimen si on alize the variables as follows: 
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where 
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Aa 2 


(5.27) 


Here, V Q is the velocity for the gas phase on the 
axis of the circular pipe. Velocity V 0 does not change 
due to the presence of the particles as the basic flow .in a 
straight pipe for the gas phase as well as for the particle 
phase is identical [^Karnis, Goldsmith and Mason (I966) ^ ♦ 

With (5.26), eqs. (5.21) to (5.25) in the non dimensional 
form become 
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(5.31) 

(5.32) 


where 
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K is 'the measure of the curvature of the pipe and it 

V a 

also varies with flow Reynolds' number R e ( = ~^*“) * ^kus K. 
is zero for a straight pipe. 

In the present investigation, for the sake of simplicity, 

a is considered as constant. Then eqs. (5.28) to (5.32) 
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Hero, eq. (5.33) Is obtained by eliminating (v-v ) 
from eqs . (5.28) end (5.31). 


If the pipe is straight i.e. K = o, then <f> = <f> =0 

P 

is a solution of eq. (5.34) and eq. (5.33), then becomes 
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(5.37) 


Eq. (5.37) is the equation for the flow of liquid 
under pressure in a straight pipe. Eq. (5.33) for v can 
be regarded as determining the deflection of a flexible 
membrane by a normal pressure, it is then clear that the 
pressure in the central part of the membrane will have the 
most important effect on the deflection. Near the center 
of the pipe, the secondary motion is roughly in the 
x-*direction, the velocity v being right handed about OZ. 
This suggests replacing the actual secondary motion by a 
uniform stream in the x-direction . I/e choose 

* = -2ky (5.38) 


where k is a constant. The boundary conditions at the 
surface of the pipe are not satisfied by the assumed 
secondary motion. 

Eq. (5.36), on using (5.38) becomes 
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Using Kongo's method., we get a solution of the 
eq_. (5 .39 ) as 

I" a 

‘L = ( + Ce Xx ) y 


(5.40) 


where G and X are arbitrary constants* 


For 4> , given by eq_. (5.4o), two cases are of particular 
P ■ 

interest, namely: 

Case I 

J>! 0 

| — j » Ce Xx 


which gives 


F* o 


<f> % --r-^ 


Thus <f>p, for X negative, is of the same nature as 4> 
i.e, secondary flow for particle phase can be replaced by 
a uniform stream. This case corresponds to the consideration 
of fine particles. 


Case II 

F' a 

| — — B. J « Ge Xx 

It gives 

4 , £ 0 e Xx y 

P 

which shows that stream function for the secondary flow for 
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The eigen functions of the homogeneous bound ar y value 
problem 


2 2 

3 v , Lt pv 3v 

— -p. + — “k - 21c = 0 

ax ay 


ex 


with boundary conditions 


( 5 .43 ) 


v = 0 at x = + 1 and y = + b 


are 


(5.44) 


cos (m + ~) , m = 0,1,2,... 

These eigen functions are useful in the construction of 
Green’s function. 


Let G(x,y;5, n) be the Green’s function for the 
homogeneous boundary value problem defined by eqs. (5.43) 
and (5.44), then G satisfies the equation 

£-! + _ 2k IS = -6(x-e) 6(y-n) (5.45) 

ax ay 

such that 

G = 0 at x = + 1 and y = + b (5.46) 

for each fixed x, we expand G in a series in eigen 
functions cos (m + |) and the coefficients then depend 
on x, S, n . To simplify the notation, we suppress the 
dependence on x,£, n 

Let G = 2 G cos (m + ■g-) ^ (5.47) 

m=0 m 


where 
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J- 

G m = b G C0S ( m + |) F d y (5.48) 

Srnce G = 0 at x = + 1 and y = + b, eq. (5.48) gives 

= 0 at x = + 1, m =0*1,2,... (5.49) 

multiplying eq. (5.45) by i cos (m + jy) ^ and then 
integrating from -b to b, we get 


a 2 & 


yr - (“ + §) 2 $ a m- Sk Sr = - k 0 ° 3 < m+ S ir ax-« 


dG. 


ax 


b 


(5.50) 


Self adjoint form of eq. (5.50) is 


2 

. ~2kx dG v b 0 , 01 

L, ril m - . e" 21oc a = 

5x L 'COs u n 3x - J cos y n m '• Ji 


oos b m n ax 


m 


(5.51) 


where 


% = 003 < m + l> § 


In order to find G- given by eq. (5.51), satisfying 
the boundary condition (5.49), again the method of Green's 
function is used. Hence the problem of constructing two 
dimensional Green's function has reduced to the construction 
of one dimensional Green's function. We construct now, the 
Green's function for the homogeneous boundary value problem 


d... r* b 
3x 


c 


cos v m n 


e ax ^ s 

e J “ ' 


ax 


cos u m ri 


m 


0 


(5.52) 


satisfying the boundary condition in (5.49). 
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ii & m ( x » ) is the required Green’s function for 

eq. (5.52) satisfying (5.49), then 

■u 9V SG b^ 2 e“ 21cx 

f- X C-o-“- m n e x .(,-,) (5.53) 


where 


= 0 at x = + 1 


The fundamental solutions of the homogeneous part of 


eq. (5.53) are 


where 


1 jr ^ ^ V 

e ^ sinh a x and e cosh o^x 


» =r?" *1 

m ’ m 


G (x, V) = e xx (A sinh cr^x + B cosh o m x) for ~1 < x < 5’ 
= e^ x ( G sinh a x + B cosh ox) for C'< x <_ 1 


(5.54) 


such that 




V 4 ) 


WD 

1 00s "m" 

G ;(«D = - =SFs T ~T 


where G^ is the derivative of with respect to x. 

Eq. (5.54), satisfying the conditions in (5.55) gives 
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n* / >- f \ COS V ™ V t V i r| ^ 

G m (x,5 ) - -rp TnTTr e U sinh o (5 r -l)sinh a (x+1) 

ill m Ui 111 

f or -1 < s < 5 1 


cos n 
in 

'BV’s rnlTTa*” 

m m 


e k(.*,+5 ) ginh a m (x~l)sinii o m ( £+1) 
for < x < 1 (5.56) 


The solution of eq. (5.51) is 


G m = { ^(^-C) e~ 2kC ' G m (x,? f ) dS» 


cos u n % 

or G m = ~ W x pZKT<T e ~ ~ X s;blh ° m (5-l)sinh ° m (x+ 1) 

for -1 < x < ? 


cos v m n 
T?3 m smlTXc m 


e “^(^“ x )silih a (x-l)sinli o m (5+l) 


for C < x < 1 


(5.57) 


Hence eqs . (5.47) and (5.57) give G, the complete 
solution of the system of eqs. (5.45) and (5.46) 


5 .31 lor circular j3r_o_ss~sect ion 

Suppose that the cross-section of the pipe is a 
circle of radius a and let x = r cos 0 and y = r sin e 
so that r = 1 is the boundary of the section. 

In eq. (5.42), we substitute 

v = ue kx (13.58) 
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and get the resulting equation as 




d~u , d~u _ „„~kx ' x 


C 2 x 2 (.— J2. + Ge Xx^ e (2x~lc)x 


, _ . -i- — .. 

sx 2 ay 


k u = -4e 


K 1 / 2 [2C 2 X e 2 Xx ~4k0 2 l!' 1 -21 > 2 a 2 /x ^ 1//2 


Changing it to p olar -coordinates, we get 

JLJi -|- A AS + A_, jLJi _ k 2 u = - 4e‘~^ :r COs ® 
ar 2 * ** r* S0 ¥ 

G 2 x 2 (~ , - < i- + Ce xrc0se ) e ( 2X - k ) rc0se 


IC 1 / 2 [|C 2x e 2XrC0se --4ka 2 l! lt - 2P 1 2 a 2 /x J 1//2 


(5.59) 


The ej^en functions of the homogeneous boundary value 
problem 


with 


are 


d 2 u , 1 5U 1 5 2 u v 2„ _ 

— o + — tz ; + “"w T ~ m u - (j 

ar 2 r 5r r 2 ae 2 


u = 0 at r = 1 


(5.60) 


( 5»61) 


cos ne , n — 0 jl» 2 , ... 

If Gr(r,e s r 0 ,e Q ) is the Green's function for the 
boundary value problem defined by eqs. (5.6o) and (5.61), 
then 

2 - « . ~ ~ 2 

L, 

6r 


sf§ + 1 a£ + 1 i| _ k 2 & = - | 6(r-r 0 )6(e^ 0 ) (5.62) 

...2 r ar QQ ^ 


such that 
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Let 


G = 0 at r = 1 


G = 2 G cos ne 


n=0 


n 


where G's are the functions of r.r and 

i i ' r) 


o 


are given by 


and 


1 ^ 

G n = - f G cos nede, n=l,2,3,.. 

~TC 


1 ^ 

G o = Sa f G cos ne<ie » n=0 

— TC 


Eq. (5.65), on making use of (5.63) gives 


& n = 0 at r = 1 for all n 

1 

For n ^ 0 t multiplying eq. (5.62) by ~ cos 

It 


then integrating from -it to it, we get 
3 Gr 2 

Ir— — + ic 2 )G = - - 6(r-r )cos ne 

6r *= ar -J v r ' n it v o' o 


The fundamental solutions of the equation 

2 


Sr O <f + fc2 > G n = • 


are the modified Bessel functions, I n (kr ) and 1 
given by 


n 


v nso^L 


m=0 2 2m+n m!(m+n) l 


and 


(5.63) 
( 5 .64) 

and 

ft 

(5.65) 

(5.66) 
ne , and 

(5.67) 

- n (L*), 

(5.68) 
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K n (kr) = (_l) n+1 l n (kr) log(jp) 


,,/kr, 


+ 


+ 


CO 

(_ i v JL.. /kr s2m p £ 1 ( m +1 ) r'^m+n+l 1 

^ 2 j 2 m=0 rnffm+nl!^ } LT(B+TJ + TTmfe+il 


1 


n-1 


2 L-ljfLMrlji (kr^m-n 

ml ^2 > 


m=0 


(5.69) 


K n (kr) is infinite at r = o, therefore I (kr) is 


n 

taken as the solution for r < r , Hence 


G n tr ’ r n) = a i(InO“)) 


n'Ti f or 0 < i < r 0 

a n I n (kr)+p ri K ri (kr) for r n < r < 1 


n n' 


(5.70) 


Conditions on G „(r,.r Q ) are 


G (1) 
n v ' 


= 0 

G f r o+> = V r o-> 

s A< r oh " G i( r oJ = - cos ne o 


(5.71) 


where * denotes the derivative with respect to r. 

Eq. (5.70), with conditions in (5.71), and using 

i 0 G„( ta o) *i0“o> - = =™ s t. = 


gives, for n 0 
cos n0 

C 


n'-~ 3 ~o / ~ Exr Tk): ^n (k ^ (kr o)"y ^ 0 ^^) 0 


G -n( r » r n) = 


for 0 < r < r 


o 


cos ne 


kitl 




for r 0 < r < 1 


(5.72) 
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Along the same lines for n = o, we get 
I (hr) 

G o (r>r o ) = ST Jk)' L I 0 ( k )K 0 (lcr o )~I o (l:r o )K o (k) J 

for 0 < r < r 
— o 

= siar o TO' C I o( k ) K 0 ( kr )~ I 0 ( kr ) K o (lc )3 

for r < r < 1 
o — 


(5.73) 


Hie complete solution G of eq, (5.6o) satisfying 
the boundary condition (5.61) is now known from eq. (5.64) 


5.4 S_o_lu__tion_ of_ the. .problem 

5 .4a ffor re ctangul ar cross-section 

Solution of the eq. (5.42) subject to the boundary 
conditions 

v = 0 at x = +1 and y = + b 


is 


1 „ 2C e x 2 (- F ll + Ge«) e 2>5 

- 1 “’ b K 1 ' /2 {S0 2 ie 2XE -4ko|i' , -2I' |2 o 2 /x[] : ^/ 2 


* G(x,y; 5, n).Vdf d n 


(5.74) 


The solution corresponding to the first term in the 
integrand in eq. (5.74) represents the case of clean gas. 

Using the identities 

cXr r - 

(i) / e a;s: sinh bxdx’= - t t — ( a sinh bx-b cosh bx) 

a ~b 



(ii) 


o p 

y~-b 


00 4(~1) cos y y 

m=o tv 3 

m 


the solution corresponding to. the clean gas 


p p v, r 00 (~1) sinh k cos u y. sinh a 3 

2( h ~y ) + 8e c — JK. _J!SL 

m=0 by sinh a 

m m 


Ti-L ' | 

(-1) cosh k cos y y cosh ax 

«, * • ^ JW .v*-.r^L ■ ■, 1 

K 00 S h % 


(5.75) 


This solution is the same as obtained by lean (1959) 
by the method of separation of variables. 

The solution corresponding to the second term in the 
integrand in eg. (5.74) represents the case of the particulate 
flow. Here, we integrate this term in the following cases* 


Case I 


pi a 

» 0e Xx , x negative 

X 


Then the expression 


2q s a 2 (- — + Ce X? ) e 2X? 
X 


2 ^ /IV 2-ra I _OV I 2 2/ “I 1/2 


X e^ -4k0^* -21* aVx “ 


9 9 .pi 0 If \ ,2X5 

2C X (—■ + — • e * ) e 
a cu 


■»: «*w. 


jT^PUSk+P’/x) 3 


TfT"' 


after neglecting 


second order terms 
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with tnis simplification, the solution corresponding 
to the particulate flow, given "by the second term in the 
integrand of eq. (5.74), is 


DO 

2 

m=0 


4-D (~l) m f'C^X cos mve 


kx 


__ _ m w _ _ __ 

hy m sinh 2a m t2KI ,t (2k+ j~(2X-k) 2 - ”J 


oo 

2 


4S (~l) m C 3 X 2 cos y m y e kx 


m= 0 ha p y m s 


ILL 

3inh 8 % [-2EI? , (Ek+ E-) ]V C(3x-kr- B ‘] 


“> 4(-l) m 0 8 XP' cos u y e 2Xx 

m “0 b^L-SKSySk+P'/X) : l7S CT2X-k)' ? - 4] 

4(-l) m 0 S X Z cos V e 3ix 

— Z ^ ( 5 . 76 ) 

m=0 bo pUm E-2Kl'(2k+P'/x)l 7 hC3x-te) - o m :i 

where 

D = e “(2*~k) s jnh o m (x-l) - e( 2X “ k ) sinh c^x+1") 

and 

E = e~( 3X “ lc )slnh ^(x-l) - e^ 3X “ k ^ sinh a m (x+l) 

The complete solution v(x,y) of the eq. (5.74) is 
now the sum of the expressions (5.75) and (5.76). 

The total rate of flow through the pipe depends on 

the integral of v talcen over the area of the cross— section . 

1 b 

Volume rate of flow = r t <say) = / / v(x,y)dxdy 
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After doing some algebraic simplifications, we get 


P 


■CO 

- E 


i ~ -y 


16b 16 „ 


°° <j rn ( cosh 2a - cosh 21c) 


D m=o sinh a m cosh a^ 


16F'LC 2 X 


m=0 by* sinh 2a [3-2KF' (2k+ ~\ x/ * £(2X~k) K ~ oT\ 


■“-g— 


ill 


111 - 


E 


m= 


1 6M0 5 X 2 _ _ 

0 l » p 4 sinh 8a B Q-2KP'(zk + £T) 


8G 2 X 


111=0 by 2 n-2KP'(2h + - r ):j 1/2 X £ (2X-3c) 2 - a 2 J 

2G sinh 3X 

+ p sr-1 (5.78) 

sopCcs^r-cr 

where 

L - k sinh 2a sinh 2A~a. m cosh 2a cosh 2x+a cosh(2A~21c) 


P 1 sinh 2X 


m 


m 


m 


and 


M = lc sinh 2 sinh 3A-a m cosh 2 a ffl cosh 3X+a m cosh(3A~21c) 


The first two terms of the eq_. (5.78) give the case of 
particle free fluid which agrees with Dean's (19 59) results. 
Case II 


P ' a_ 


12. « Ce Xx , 0 negative 


Then the expression 
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2G 2 X 2 (- — + ce xx ) e 2xx 


F'o 


[20 2 X e 2Xx - 4kg 2 I' , -2F t2 a 2 /x l 1 /^" 
P p x - J 


rs 


2X° n XX / ^ , °n „ W : x 3 F'a^ 2kD’ r o c 

If* ue ( "T ' + Je ) SK 0 (’^4 + 'X™*') 


• 2 2k? r " 2 

> 

r + ”■5’ 

CTx" CTx 


I"a 

( —*>‘--2. e x + 0 ) 


ii/ith this simplification, the solution corresponding 
to the particulate flow, given by the second term in the 
integrand of eq, (5.74), is 


~ 2(~l) m f 2 X^/k ce kx cos P rfl y F'o 

m=° bXv m [;(X-k)^- ^fjsinh 2a m 

x £e“( X '" k ) sinh cr m (x*-l)-e^ X “^sinh a m (x+l) 3 

“ it -r- ®1° £ 003 v e xx 
m=o bxy n [;‘(x-k) 8 -%ii 

~ 2(-l) m f 2X 3 /EC e kx 0 2 cos v m y 

m =0 *% £ ( 2X -k) T - sinh 2c m 

x [V (2Uk )sirh a m (x~l)”e^ 2X ’" lc ^sinh ° m (x+l) J 

- sx-iFu* s /& ° 2 005 y e 2 x x 

m= 0 i>y m L (2X-k) 2 - crJ5 

2 2 (- 1 ) m ^ 5 / 2K Q 005 V_ll. p , a e - xx 

m=0 hXy m £(x+k) 2 “ P 
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“ i ( ;. 1 . ) _ m .. 3? .'!£. T j7 2K °£! fi_ 00s v 

m=0 1 ° x \ l L (X +lc)^- 0 2 ~] sinh 2 ^ r 


kx 


m 


x j-- e ( x +k) sixth o m (x-l) - e“^ +k ^ sinh 0 m (x+l) 2 


CO 

+ s 
m=0 


2(-l)“ 0 2 lfx T /2Ii e kx E 1 OOS u y 
s “ h 2o m 


x (e lc sinh a (x- 1 ) -e~ k sinh a (x+ 1 )) 


+ 2 
m=0 


00 2(~l) m C 3 V" X^/2K cos P m yP 1 


bn 


T" 


m 


(5.79) 


where 


P«a 2 Ska 2 !' 1 
p , p 

— —ij—fc + “~w— “ 
1 0 Ki x 


The complete solution v(x,y) of the equation (5.74) 
is now the sum of the espressions (5.75) and (5.79). 


In this case, 

• 1 b 

Volume rate of flow = Pp(say) = J f v(x,y)dxdy 

^ -1 ~b 


p = 16b 3 _ 16 ~ *4 00ah 

3 m=0 yf sinh o m cosh c m 

oo 8 f 2X 3 /K P ! 0 0 
m=° by^xQx-K) 2 - a 2 jsinh 2 

X U cosh(x-2k) + k sinh 2 <y 3 inhx~ 0 m cosh 2o m cosh x“] 
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CO 

- 2 


Si f 2A 3 /K G 2 


m=0 by* A[~(2A-k) 2 -a;H gfcfh 2o' 


m 


oo 4 fSX^/K IM c 

m=0 bAp^ [ (A+k) 2 - sinh 2a 

in • m 


m 


oo 

+ 2 


* !-■ a iu cosil ( x+2k )“ k sin51 2a m sinh A-a m cosh 2a m cosh A^ 
4 lf2A 3 /K C E |(a m cosh 2k-o m cosh 2o m ) 


m=0 

bM m sinh 2o m 

OO 

L! V2 X^/K 0 sinh X 

ju 

m=0 


oo 

8 f2* 3 /K G 2 sinh 2 A < 

2 

m=0 

2hXM^ C(2X-lc) S _^J + m 

00 

y 

8 f A 3 /2K f « a C s inh * 

u 

m=o 

hx“ £ C(x+k) g - °IJ 


fey 


m 


( 5. 3q) 


The first two terms are same as in case I, giving the 
volume rate of flow for the particle free fluid, same as 
obtained by Dean (1959)* 

5 ,4b p^oular__c r os s -s e ct i on 

The solution of the equation (5.59) is given by 
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1 n kr cose 
u(r ? Q) = J J {4e 0 0 

0 -71 


G(*,e;r 0 .e 0 ) r 0 to o ae 0 


i n . c 2 x 8 (I'> + Ce ir o 00se ° )e (2x - lc)r o oose o 

0 ' IcV^tJO^e 2 ^ 000800 - 4ka 2 

P P' - 1 

G(x,ej JC 0 ,e 0 ) 21 r 0 dr 0 d 6 0 (5.81) 

Tae first term in the integrand in the eq_. (5.81) 
gives the solution for the particle free fluid, 
using the relations 


% -kr cose 

/ e cos ne o de o = 2 ti(- 1) I n (kr Q ) for n > 0 

—71 

(5.82) 


/ rl^kr^OcrJclr = f-jilykr) I»(kr) - (kr) I^(kr) 


+ Ek V 1 ®) ^>3 


(5.83) 


/ rl~(kr)dr = g~ Cljkr) I, ^( kr ) - (I^(kr)r 

+ |e Ii(kr) I n (k) 2 (5.84) 

I he solution corresponding to the clean gas flow is 


co (_l) n cos ne I* (kr) I (kr) 


21 (kr) I’(k) _ \ — ~n' 

.0 + 4 2 

kI o (£ ) n=l -n 


2rl l (kr) (-l) n r cos ne I'(kr) 

r:— | — + 4 2 t ^ 

n=l 


(5.85) 



The second term in the integrand in the eq. (5.81) 
gives solution corresponding to the particulate phase. 
Here, we will solve it for the two cases. 

Case, I 


l"o 


» Q e ^ rc °s0 ^ ^ negative . 
X 


The expression 


2 9 p ' a n ^ r n COS0 n (2x-k)r OOS0_ 

+ Ce 0 °) e 0 0 

A 

[$C 2 xe ^ 0 S9 ° - 4ka 2 I ! ~ 2F ,2 c 2 /x[] 1 / 2 

P P 


pr (2x -k)r Q cose o 0 (3X^k)r o cose 


f e 

!el. 


o 


-2KF ’ (21c + P'/ X )J 1//2 


Then the solution of the second term in the integrand 
of the eq. (5.81), after using the eqs. (5.82) to (5.84) is 


o 2 xi'> £i 0 (ic)i 0 ((k~ 2 X)r)-i 0 (kr) Ip( k ~ 2X )H 
I Q ( k) [ 7 2KF ‘ ( 2k+F ‘ /x ) jk^-(k~2 x) 2 fj 


+ 


=o SC 2 XP' (-lfoos ne Cl n (k)I p ((lc-2 x^j ^kr^tfc-gx) J 

n =i Jjw E^s'Va+V ' A ) *? 1 



J-C> jL 


+ 2 


“ 20 3 x 2 (-i) n oos ne[;i n (k) I n ((lc-3x)r) - I (k-3x) I (kr) H 


n=l a p I n (k) I^Kf<(8k«>/») 

ihe complete solution u(r,e) of 'the eq. (5.81) 
is the sum of the expressions (5.85) and (5.86). 

It is known that 


(5.86) 


e ~kr cose = x Q ( kr .) + 2 2 ( -l) n I^kr) cos ne 


n=l 


n' 


(5.87) 


Differentiating (5.87) with respect to k, we get 


- r cose = rl^(kr) e lcr cos0 + 2r I (-l) n I^(kr)cos ne e 

(5.88) 

Using (5.87) and (5.88), eq. (5.58), with u(r,e) 
given hy the sum of the expressions (5.85) and (5.86), gives 


kr cose 


o O jfcrcose 03 

v = — 2 a I (kr) cos ne 

k K n=0 


o° 

0 2 xp » gkreose £ e (2X-k)roose„ £ p n I n (kr) cos ne J 


+ *~-*--* 


[-210?' (gk. + r/x) 3 A/<S [k -(k- 2 *) I] 


n=Q 

.2 / -i n i \ 2' ” 


n 3 x 2 e kr cose r e ( 3x-k)rcose_ 2 f I (kr) cos ne J 
u n=p n 

+ “T" , fS5’«*(2k + i'/x) :] 1//& L k2 -( k “ 3x ) 2 5 


( 5 .89 ) 


where 


a 


0 = ij,(k)/i 0 (k) . * n = 2(-if iyk)/yk) 
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I 0 (k-gx)/l 0 (lc) , p n = 2(_l) n I n (k-8X)/I n ( !c ) 


= I 0 ( k - 3X )/ I 0 ( ,!: ) > \ = St- 1 ) 11 I n (i- 3X )/I n (k) 

1 It 

Volume rate of flow fg(say) = / / v(r,e) rdrde 

0 ~7t 

Using the identities (5.82) to (5.84) and 
1 I-i (k) 

/ rl (kr)dr = - ir ~ (5.90) 

o J - 


we get 


v'w 

l { I?, a. 


TtO^XP' p. 


n 


n=0 


n 


+ i"/ x )H 1/,2 Lk z -(ic-2x} 2 3 


itC 3 x 2 r. 


n 


a {t— SICI? » ( 2k + 3? ' • /x ) Q?~(k~3 X) 2 J 

P 


* £ (l+|“g') i 2 (h) - n 

7iC 2 X]?’I 1 (2X) 

h X JU-2KF 1 ( 21o+P 1 / X ) □ 1/2 nin 2 -(ln~2 X)^ J 
2tc0 3 X 2 1^(3 X) 

+ 3X a p [-20 1 (2M'/x) H r/ " 2 Lh 2 -(h-3 x) 2 5 

(5.91) 

The first term of the BUS of the eq. (5.91) gives 
the rate of flow for the particle free fluid, which agrees 
with the analytical results of Dean (19 59) » 
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Case II 


I ' 1 a 


-£■ « oe Arc ° s 0 


) C negative 


The expression 


A 2 + Ce Xr ° QOSG o- 1 p ( 2x " k ) ^ 0 cose 

-TaFoo^ 


£2C^Ae 


> “'4kcr^E , .»- 2P> 2 a 2 /A“] 1//2 


C 2 a 2 (-— • + oe xrcose ') e xroos 0 i CA 2 !! - D 

x ^7?7S ( ~ 


P 1 a 


p» a 


■x e 


-2Ar cose 


V/ith eqs. (5-81) to (5.84), eq r (5.58) give; 


v 


2r cose , 2e kr cos ® 00 




+ 


T 


2 0 n I n^ kr ) 00 s n0 


n=o 


4 - 


0 2 A 2 P«c e kr COse 
p 

[Ska) 1 / 2 A^l(k-x) 2 -j 


[e(A-k)rcose_ 2 a^I 


n=0 


q 3 x 2 e kr cose 


4 “ 


— sr-m V* — — ( 2 x -k)rcos0_ V t) '] 

( 2K X) 1 / 2 (k 2 -( 2 A~k) 2 2 L n=0 n 


C API' 1 o e kr COs 9 1 , 

+ — j^(2x+k)rcose^ ^ c 

2(2Ka) /2 [Ic -(2A+k) 2 J n=0 1 


C 


2 x 2^ e kr cose 


+ jx( A+lc)rcose y d j 

2(2KX) 1 / 2 f T -(A + k) S 3 n=0^ r 


0 


Ce» r oose ) 


; n (kr)cos nej 

; (kr)cos hqH 
n v 


I n (kr) cosne^ 


(kr)cos ne_ 


where 
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V s I 0 C»-J c )/ I 0 (fc)i a a =SC-l) n I n (x^:)/I n (k) 

V l 0 (2x-l:)/I 0 (k),b n =8(-l) n I n (2A-i)/I n (k) 
c o“ y 2 xik)/I 0 (k),c n =2(-l)\(2Art)/yk) 


V I 0 (x+k)/I 0 (k), 


1 n 

Volume rate of flow = !4(say) = f f v(r,e) rdrde 

o -it 


,2. 


F 


.En 


*C XF' c a 

a — 1 ( £r- a - — T7r^Tr - ' 

*** ic n / pit • 


(2KX) 1/2 6 i: -(X-*) 2 n 


lt C 3 X 2 b 


*°¥' IJjL 


( 2XX) 1/2 2 ” 2(2KX)^§ Z -(2X*fcj S D 


— rq > C(i+ « 3 

2(2IC X) T/2 [jp-( X+kf n k 


+ 


2%0 2 I'o I x (X) 


itCTX I 1 (2X) 


( 2K~Xp'^^-lTI73 + (2KX) 1 ^: 2 ^(2X~lcl r 3 


HC^P , o p I 1 (2^+ 2b ) 

+ rr^Tfer^' — 


;_T2 -i 


«^(H k ) &•<**)* 1 

tiC 2 X 2 P i I 1 ( X+21c) 

(2K X)^(x+21c) ^Tx+fc)^ I 



J -OU 


5 , 5 hi .a cu s a i on,, of, res ults 

In this section, we present the certain basic features 

of the gcs-p articulate flo w through a curved pipe. . I' or 
numerical computation, the following values of the parameters 

are considered: 


k — 0 , 1 , 2 , 3 , *....,10 

c = -.25, .25 
X — — .9 , »9 

3?» = .05, 10, 15 
K - 1.0 


°p = * 6 » 1,5 

pig. 2 gives the variation of the volume flow rate with 
respect to the secondary flow parameter k, when the secondary 
motion for the particle phase is a uniform stream, fig. 2 
indicates that as k increases, the volume flow rate decreases 
Also, the volume flow rate for the clean gas is more as 
compared to the flow rate for gas particulate system. This 
is due to the fact that relaxation time is small, so that 


the particles are fine. 

Pig. 3 gives the variation of the volume flow rate with 
the secondary flow when the stream lines for the secondary flow 
of the particle phase are approximated hy the exponential 


jurves. fig. 3 


in this case also, volume 


flow rate d< 
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clean gas, volume flow rate is less as compared to the flow 
rate for gas-particulate system. It is due to the fact 
that relaxation time in this case is large, so that the 
particles are coarse enough. 

Pig* 4 gives the variation of velocity v with the 
secondary flow at y = 0 for case I. This figure shows 
that v reduces as the secondary flow through the pipe 
increases . This figure further indicates that the region 
where the primary motion is maximum, shifts outwards. 

Also, the velocity v for the clean gas is more as compared 
to the gas-particulate flow . 

Pig. 5 gives the variation of velocity v with the 
secondary flow at y = 0*0 for case II. This figure indicates 
that in this case also, v reduces with k and the region, 
where the primary motion is maximum, shifts outward* Further, 
the velocity v for gas— particulate flow is more as compared 
to v for the clean flow. 

Pigs. 6 and 7 give the variation of volume flow rate 
and velocity v with the secondary flow parameter k for the 
oase I. Por different set of parameters, same trend is observed 


as in Pigs • 2 and 3 
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CHAPTER VI 


■'*- ■ ***' «r. *• 4fcra»va«?,r. art- /deutsfla. aero* 

GLOBAL STABILITY AND UNIQUENESS .OF AIT INITIAL 

BOOT NARY VALUE PR03LHI OP A LUSTY GAS MODEL 

% 

6.1' Introjiuction 

Two types of uniqueness theorem of an Initial boundary 
value problem (IBVP) for Navier-S tokes equations have been 
pointed out. First of them ensures only one solution for 
each initial data while the second guarantees only one steady 
state solution for not too small a viscosity. These results 
are of significance in hydrodynamic theory of stability. 

F oa (1929) found first type of uniqueness theorem for 
initial value problem. This theorem ensure , uniqueness of 
regular solution of initial value problem, provided such 
solutions exist. If solutions develop discontinuities, 

Foa* s theorem may not hold. It is generally supposed that 
discontinuities can not develop in a viscous fluid but 
despite great efforts by eminent Mathematicians, It has 
never been proved [Ladyzhenskaya (1975) yj , Serrin (1959) 
while discussing the concept of universal/Global stability, 
pointed out the uniqueness of steady solutions for not too 
small a viscosity. In fact, the two types of uniqueness 
theorems result as a byproduct of Serrin* s concept and 
analysis of universal stability [Joseph (1976) [[] • 
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Serrin (1959, 1963) extended the uniqueness theorem 
to the case of compressible fluids. However, universal 
stability theorem is not available till today in such a 
s ituat ion . 

The purpose of the present chapter is to search for a 
uniqueness theorem and a theorem for universal stability for 
a dusty gas. model. 

There are certain works in this direction. Crooke (19 76) 
examined a uniqueness cr iter dm and Dandapat and Gupta (1976) 
obtained a universal stability theorem for a dusty gas model. 
In both the works, the authors considered the constant number 
density of the particle phase, Dandapat and Gupta (19 76) 
did not obtain any uniqueness criterion while Crooke (19 76) 
did not bother about universal stability theorem. The two 
works are quite independent. 

In the present chapter, we improve the universal 
stability theorem of Dandapat and Gupta (19 76) and as a 
consequence, we establish a uniqueness criterion. At one 
point, we depart from the analysis of Dandapat and Gupta to 
obtain a much improved theorem on universal stability. It 
ensures a wider range of stability region. ¥e further 
generalize the analysis to cover the case of non-uniform 
number density of particle phase. One is able to extract 
only a uniqueness theorem. 
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6 *2 Fo rmula tion of the problem 

Consider a viscous incompressible two phase flow in a 
closed bounded volume ' y(t) whose boundary is ?(t). We 
assume the gas phase to be of constant density while the 
particle phase to be incompressible with non-uniform number 
density. The fluid is set in motion by external forces or 
by the motion of the boundaries. If U(X,t), V(X,t) defines 
the velocity fields for gas and particle phases respectively, 
x(X,t) the pressure in Y(t) and K(X,t) the number density 
of the particles, then basic equations are 

gi + (U« v )2 “ + vx - F(X,t) - Sjjr (Y-U) = 0 (6.1) 


div TJ = 0 (6.2) 

SY y 

gi + (I* V )I = s (2-v) + 1 (X,t) (6.3) 

|| + (V. v) H = 0 ' . (6.4) 

div Y =0 (6 .5) 


Here, X S Y(-fc), t > 0. F(X,t) denotes external 
force per unit mass, v the kinematic viscosity, p q x the 
pressure, m the mass of each particle, If the drag coefficient 
on the particles. We assume the basic equations to be one 
as developed by Saffman (1962). 

Further, we prescribe the initial and boundary 


conditions as follows: 
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u(x,t) = u s (x,t) X e 

5(t) t > 0 

(6.6) 

*1(2, t) = n.Y s (X,t) 

xe ? (t) t > o 

(6 ;7) 

U(x, 0 ) = U 0 (X) x e 

Y ( 0 ) 

(6.8) 

V(x,0) = V 0 (X) xe 

Y(0) 

(6.9) 


where 

div U = div V =0 
o o 

y(t) denotes the volume occupied by the fluid at any 
instant and n the outward drawn unit normal at a point of 
C(t). 

Since the two phases are assumed to he incompressible, 
total volume remains invariant, however its shape can change . 
fhus 


d 

dt 


/ d T 
Y(t) 


where we define 


= J div T(X,t)dt = 0 

Y(t) 


d_ = 
crt Qt 


+ (Y.V) 


It is to be noted that material surface is defined with 
reference to particle phase. Equations (6.1) to (6.9) 


define an initial boundary value problem (IBVP). 


let U(X,t,U ) and V(X,t,V 0 ) denote a solution of 
(6.1) to (6.9) for a given data { U g (X, t) , Y g (X, t), l(X,t), 
Y(t)> for a fixed v , K,m, p . The main problem is to 
establish two theorems of the following type : 
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* ^ or e ach U 0 , V Q , there is only one XJ(X, jU Q ) 

321(1 satisfying IBVP (6.1) to (6.9). 

The or era _ 2 . If the given data { U , V ; E, y } is steady and 
if v > v then there is atmost one steady solution of I3VP 
( 6 .1) to ( 6 .9 ) . 

As a byproduct of the stated theorems, we will obtain 
a global stability criterion. theorem 2 is essentially 
significant in the sense that it ensures one steady solution 
provided kinematic viscosity is large enough. For a smaller 
viscosity, problem of turbulence arises and for a larger 
viscosity, fluid may not remain Newtonian and henoe basic 
equation (6.1) may not remain valid. 

To prove these two theorems we build relevant structure 
in the following sections. 

6 .3 Ene rgy equ ation 

let (U X,t,U 0 (X) > , V' {X,t,7 0 (I) > and U* {X, t,U 0 +u Q (X) >, 

IP (X.t.V + v (X)> be two solutions for fixed v «K,m, pand 

same data {U (X,t), V 0 (X,t), E(X,t), y(t) > . Denote 
"“ J S *■"* ““'S *“* 


ii(X,t) 

= U’(x,t) 

- E(x,t) 

l(X,t) 

= r(s»t) 

- v(x, t) 

P (X, b) 

= ix 1 (X, t) 

- 7i(X,t) 

n(X,t) 

= H'(X,t) 

- N(X,t) 


Then. the bo.sic equations (6.1) to (6.9) reduce to the following: 
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3U 2 

+ (U.Vu) + U, V U + U.yU “ V V B 



+ VP - (v-u) ~ ( v-u) - ~(7 -IT) = 0 

P P 

(6.10) 


div u = o 

(6.11) 

82 

St 

■ ' K 

+ Z • V Z + 2 • v 2 + 2 • V X = ~ ( B“X) 

(6.12) 


div v = o 

(6 .13) 

8n 

sT 

+ (v*v)n + (v.y)l + (f. 7 )n = 0 

(6.14) 


u(X,t) - n.v(X*t) = o ie C(t), t > 0 

(6 .15) 


a(2,0) = u Q (X) X 6 y(0) 

(6.16) 


2 ( 2 , 0 ) = v 0 (X) 2 e y( 0 ) 

(6.17) 


Define 


E(t) = §< |u| 8 > = | y / t) |u|V 
e(t) = | < |v| 2 > = | y [ t ]2l 2d ' 


It is to be noted that 


It E ^) = I v ( t) Ie < hl 2 > dI + I y ( t) ll 2 ^ 


1 r d 
8 T (t) ^ 


av. 


Hl 2 d -| Y { t) bfe 3 lHl 2 d -l Y { t) h. ^ 


.1— ( Y . j v 

P. ir V i I — 


2 )dx 


1 f ft' inl 2dT + | / .gtr* 

2 Y (t) ot 2 Y(t)°-3 

r / k llil 2 at +| / n .v .|u| 8 as 

2 y( t ) at 2 c(t) 3 3 

- S u fl dT * 

Y (t) - Qt 


dr 
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The last integral vanishes owing to first condition of (6.15). 
further 


e(t >=5 L, I? lh 2dT + % / Jv| 2 aiv 7 dT 


Y(t) 


y(t) 


Thus 


dv 

/ v « d? 
Y ( t ) ' dt 


= / v (“p v + V.vv)d- 

Y(t) ~ dt “ “ “ 


KJ 

E( t) = f u# d T 

m T(t ) -■ ° l 


(6 .18) 


&’<*> ■ Y ( t) Z-<fc* + S' v S> dt 


( 6 ,19 ) 


Making use of eqs . (6.10) and (6.11), one therefore obtains 


4r E(t) = / [ 

dt Y(t) 


s 2 u, ?n au. 

vu. — — u. r^r" — u- u. vtr~ 
■ l 8 2 l 8^ l 3 ax 3 

J 


5TJ. du^ 

~ u i U 3 ST! “ U i U j 9Xj* 


+ (u 1 v.-iu| 2 ) + | n( UiV^jul 2 ) 


+ V ("iWi> 


Y { t) alr ^ 


au. au. . 

_ r r i)dt - / dV.u u dT 

y( t ) ^ X j 8X 3 v(t) 13 1 3 
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where 


+ s / (H+n)(u T.-|u| 2 )a T + I / n(u 1 V u i U i) d ' 

Y(t) Y(tJ 


, 5TJ . dU, 

hj = s ( bit 3 + 

Use of boundary conditions on v and Gauss theorem 


yields 


E(t) 


■ 5U 1n2 


- / P D. .U.U.+ v( v --i) 3 d -' 

Y(t) U 13 1 3 

+ £ / n (u.V.-u tipat 
P 1 X XX 


+ I / (U+n)(u.v ~|u| 2 dT 

Y(t) 


( 6 . 20 ) 


If particle density H -happens to be a constant so 
that n = 0, it reduces to 


5 u i g 

E (t) = - _/ cuuiv v W 34t 


Y(t) 


+ ~r / (Up^-iuj )d' 

Xt) 


(6.21) 


Similarly, if use is made of eqs. (6.12) and (6.13) 
then eq. (6 .19) gives 

4 e( t) = / v. E-t.vT- 2 -vy + | (£-2) 3^ 

h t: x 7 /i\ 


y(i) 


= - / ( v j^v . f~ : l t + Va ^ 

Y ( t ) J '* 3 3 


+ ! / (u i v i -|vj 2 )dT 

m y( - t) i u 


v-)dT 
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In this section, we state/obtain few subsidiary results. 

A 

Lem ma 1 . If L m < 0 is the smallest of the three eigen- 
values of jjl.. J over y( t ) , then 

- <u ,D «u> = | f> m j < u j 2 > 

Proof . I (2,t‘), o < t 1 < t is a symmetric matrix defined 
over Y(t'). Por each fixed t ! and at a point X of y(t’), 
let D 11 (X,t>), D 22 (X,t‘), D 33 (X,t‘) be the real eigen values 
of fKlj-tl). But since 


T r S = div = %1 + J 22 + J 33 = 0 
therefore, one of the eigenvalues is necessarily negative . 

Let < 0 denote the smallest of the three eigen 
values over y(t') up to time t. Then 

An = ^11 ,:D 22» :D 33^ 

0<t'<t 
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low at a fixed t' and a fixed point X of Y(t r ), 
u.D.u = u «M D , m .u 

wr# ««., g»*J Q_ mm* mr** 

where Ij ^ is diagonal matrix and M is an orthogonal matrix. 
Thus 

u . £ . u = + D 22 V Z 2 + B 33 V® 

> III 2 min (B 11( D g2 , I) 33 ) > Bju| Z 

where 

V = M 1 . u 

and 

ill 2 = i«i 2 

Thus 

- <u. D . u > < |hj < |u| 2 > (6.24) 

This completes the proof of lemma 1. 

Similarly if d fll <0 is the smallest of the three 
eigen values of [d. . J] over y(t) up to time t, then 

- <v . d . v > < |d m J < |v| 2 > (6.25) 

It follows from the fact that div. V = 0 

Lemma 2 . (Poincare inequality) If u(X,t) is any smooth 
function such that u(X,t) = 0:X e ^("t) and div u = 0, 
there exists A > 2 such that 

4 < M 2 > .> <ii 2 > 


( 6 . 26 ) 
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where is the smallest distance "between two parallel 

plane® which just contain y(t) and 1 = mar ^(t 1 ). 

0<t’<t 

One may he referred to Joseph (1976) for the proof of 
lemma 2 . 

Use of lemma 1 and lemma 2 reduces energy equations 
(6.20) to (6.22) to differential inequalities: 

ft E (t) < 2 ( l|l ~ -p E(t) 

+ * y(t) nu i (T i^i )dT 

+ — I H(u. v,-|u| 2 )dt (6.27) 

P y(t) 1 1 " 

lor a constant number density 


:(t) < 2 ( i\l - tz> E(t) + if / (u^-iui )a 


IW 6 . 28 ) 


y( t) 


It- e(t) < 2|||e(t) + 1 Y ( t) (u l V i-hi 2)dT (6 ; 29 > 

fe shall quite often he using estimates of the integrals 


of the type 


/ F(x) g,(x)dt , F(x) > 0 over y(t) 

y(t) 

T ive can estimate it in two V'/ays • 
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.bet 


Then 


„ max It,/ , , 

* = Y(t) l F ( X )l 

o<t<t« 


2 f i (x) M x ) dt < f *<*) (f? + g?) 


Y(t) 


Y(t) 




< a / (fj + g 2) dT 

Y(t) 1 1 

Thus 

/ J(x) f^x) g i (x)dT < a (|- / f^dT + l f g?dr) 
Y(t) ii ^ Y (t) 1 4 Y (t) 1 

(6.30) 

Further, using Schwarz inequality, one gets 


( / yw f ± (i) < / r s (x)f?(x)dt / gb x )ax 

Y(t) 1 1 Y(*c) 1 Y(t) 1 


< Ct 2 / f 2 ( X )dT / g?( X )dT 
Y(t) 1 y( t ) 1 


Thus 


S ^(x) f ± (x) g ± (x)dx < 2afF/"f?dTY fjgf <1 t (6.31) 

y( t ) x x ^ x 6 i 

-p2 


< a(~ / f Jdt + | / gj dt) (6 .32) 


In fact, it is obvious from above that (6 .30 ) follows 
from (6.31) , Moreover, we note that l'(x) > 0, is not relevant 
if (6.32) the same as (6.30)3 is derived from (6.31). 


The two estimates are important and we use them in the 
analysis to follow. We list below certain estimates which 
are useful in further analysis. These can be derived along 
the lines indicated above. 
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where 


u i v i dT < e ( t ) + E(t) 


V(t) 

/ 

Y(t) 


/ ^ u.Y.dT < 2 fhjrj fifty 


f nv ® d x < p (e+e) 

v /4-\ 02. . — r l v ' 


Y(t) 


D ' 


/ nu.(V' - Ul)dT < p (e+E) 
Y(t) 1 1 < 

/ lu.v.dx < p (e+E) 
y( t) 1 1 - ^ 


9=|/ n 2 a* 

2 Y(t) 


> Pi 


p 0 = max I V* -LT* I, 

2 V(t) 

o<t<t» 

p 4 = min N 
4 Y(t) 

0<t<t 1 


max 

Y(t) 

0<t<t f 


max 

Y(t) 

0<t<t f 


fix’. 


I 


(6 .33) 
(6 .34) 

(6 .35) 

(6 .36 ) 

(6 .37) 


(6.38) 


6 .5 Mam theor em s. 

The c ase of variable numb er d e n sity 

The relevant equations/inequalities for the present 
purpose will he (6 .23), (6.27) and (6.29). The use of the 
estimates (6.33) to (6.37) for various integrals yield 


de 

dt‘ 


< .p 1 (e ' + e) 


(6 .39) 


dE / r o i-A I 2Av k fl , K P 2K R -i™ 

3E < L 2 l D J ~ fT + “ ^2 + p ^3" T 


+ 1 Pg 0 + I Ps e 


(6 .40) 



151 



(6.41) 


(6 .42) 


It is trivial to observe that {3 > o * 

Theorem 1 . Integrals of (6.42) over [o,tQ gives 

. e(t» ) + E(t* ) + e(t* ) < [’e(o) + E(o) + e(o) (pt T ) 

where E(o) and e(o) are given by 

E(0) =| / |u | 2 dt 

Y(0) 

e (°) = I lx 0 l 2dT 

* Y(0) 

and 

6 ( 0 ) =0 

In case v Q = 0 = u Q , e(t’), E(t«), e(t») will 
separately vanish since e(o) and E(o) vanish. It 
implies that only one solution will start from each initial 


Clcl/trSl • 
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T he ,, ea se of _ cons t a nt^ number de nsity 


Eq. (6.23) will "be redundant and only differential 
inequalities (6.28) and (6.29) will be relevant. Use of 
estimates (6.33) reduces the inequalities to the form 


dE 

eft 


< ( 2 A 


III 


i< 


2Av 


ffl .) E(t) •!• s e 


|) e(t) + | B 


(6 .43) 
(6 .44) 


It will he noted later that this estimate gives a 
wider range for stability region as compared to Landapat 
and Gupta (19 76) who used the estimate as provided by (6.43). 

Inequalities (6.43) and (6.44) yield 

(E+e) < p 0 (e+E) (6.45) 

where 


"s-fv + l . + (6 - 46) 


p 0 = max L?l\ 

and is a constant* Inequality (6.45) holds for all times 
for which £ is a constant. Proof of Theorem 1 follows 
exactly along the same lines as for the case of variable 
number density. 


p roo f of Theo rem 2 . Let j3 Q = -£ 821 ^ F > 0 , then 

inequality (6.45) gives 

E(t') + e(t’) < £E(0) + e (0) U e ^P (~F^ ) 0<t*<t 
Eor a steady solution E(t') = E(o) and e(t’) = e(o) • It 
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is possible iff E(t«) = e(t») = o. Thus only one steady 
solution is possible* 

One fails to prove this result if happens to be 

positive. We note that if both the expressions in the 
paranthesis of (6*46) are negative, p* will definitely be 
negative and under 'the condition only one steady solution 
will be possible. Thus 



the system. Alternatively, one might put that 

<*<£-?> < 8 (rs-l\!) ( 6 .47) 

JJ 

is a sufficient condition only for one steady solution of 
the system. It is implicit that the condition is meaningful 


provided 


and 


H . 1 
p x m 





(6.48) 
(6 .49) 


These conditions can also be regarded as sufficient 


conditions for stability for the solution U,V. Inequality 


(6.48) means that mass density of the dust must be less than 
the density of the gas and (6.49) implies that Kinematic 
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viscosity must not be too small. Thus for not too small a 
Kinematic viscosity of gas, having light suspended particles, 
(6.47) ensures stability, 

6 .6 Improved results 


We note that character of {3 in the inequality (6.42) 
and that of in the inequality (6.45) is of basic 

importance, Por uniqueness theorem, it is immaterial 
whether (3 and (3 are positive or negative. However, 
for stability theorem, their negative character is a must. 


Keeping the observations in view, we improve the 
sufficient conditions as given by (6.47) for the stability 
of the flow. We attempt to show the existence of x^ > 0 
and Xg < 0 such that (6.43) and (6.44) yield 


Ih ( e + x i E ) 5. x 2 ( e + X 1 E ) (6.5o) 

lor such x-j_ and the following must be identically 


true . 


and 


x 2 ~ 2 l d ml 


£ + x, a 

m Ip 


X 1 X 2 



(6.51) 
(6 .52) 


where 



The two equations give that 
quadratic equation 


X 2 


must be a root of the 
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nidi + 2i» m i - I-t-^3 


V 

i 

V 


+ ( g l a ml "l> C8IJ? ! - 




m 1 


KJ[ 

p 


v x A 
"!*/ “ 


v 


pm 


= 0 


Since Aq is required to be positive and Xg < 0 

One of the sufficient conditions for the existence of 

K 


Xq > 0 and Xg < 0 will be that S { cl m j ~ ~ < 0 and a 
root of the above equation lies between 0 and 2jd j 


K 

m 


Therefore a sufficient condition of stability of the 
system is that 

K 


8 KI - m < 0 


and 


( 2 idl "I) ( 2 l\i - 


m v . A 

p “ v* ' " pm 


and 


c*i«j-|) 8 -c*iaj midi + idi - 


o 


+ (2 Id I- -)(2 liL l~ ~ - -tf) are of 
pm v 1 m 1 m' v 1 m 1 p v*' 

opposite signs. 

The last of the above egressions is equal to 


EN 


A 


which is negative. Therefore a sufficient condition for 
stability is that 

2 iy-l<° 

and ( 2 l d mi - !)(2l £ - ' 

It implies that the flow is stable if 


-) _ sis. > o 

m> p v* J pm 


2 d 


m. 1 


K 

m 


< 0 
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and 


we note that the condition fails if 

we have thus established the following theorem. 


Theorem . A sufficient condition for stability of a dusty 
gas flow is that 


2|aJ-§ <0 ( 6 - 55 ) 

2|Bj ~~ " V < 0 (6.54) 


and 



A question arises: Could the analysis developed in this 
section lead to a sufficiency theorem for the case of a 
variable number density? let if possible, there exists 
G i > °, Cg > ° and Cg < 0 in. such a way that the inequalities 
(6.39) to (6.41) yield 


-| t - (E + C^e + C 2 e) < C 3 (E + C 1 e + C 2 e) (6.56) 


This differential inequality gives 


°3 “ 


C 1 C 3 


KP, 


K 


I- + % ^2 + ^3~ 2 M + m °1 

L 


m 


+ -i> + °2h 


C 2°3 = p P 2 + °2 P 1‘ 
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The last of the equations imply that CgtC^-B^ = 7 P 2 * 
lor Gg to he positive, one requires that p 1 and pg 

being positive will consequently be positive. It shows 
that Cg > 0 and 0^ < 0 are not consistent , Present 
analysis therefore does not yield a stability theorem. 


6 . 7 A comparison o f stabilit y theorems 

lor a comparative study of stability theorem, it is 
desirable to re derive Dan dap at and Gupta’s stability 
theorem . 

let us denote E = and 0 = Kg . Then use of the 
estimate as given by (6.34) for the integral / UjV. dx reduce 
the inequalities (6.28) and (6.29) to the forms 




Av 

F 


at) % + a Kg 

p ' 1 p 2 


and 


Iff K 2 i t K 1 + ( i^ml " IC S 

The two inequalities imply that 

Iff ( K l + K 2^ i Y ^ K 1 + K 2^ 


where 


Y = max [„ | D i 


m 


Av KN ^ K 

72 ~ T + m ’ 

JLi 


system is therefore stable provided 
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S + KH'i 

m p 


a + |<o iap 


£ + a < 0 
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hg_ ore m . A sufficient condition for the stability of two 
phase flow for a constant number density case is that 

• laj < K(|-r> < |£ m l) (6.57) 

jj 

It is also desirable to note the following theorem, 

Iheorem . A sufficient condition for the stability of single 
phase fluid flow is that 

~Z - l\l > 0 (6 .58 ) 

JJ 

We will now compare the sufficiency theorems for 
stability of two phase flo.w as provided by the inequalities 
(6 .47), (6 .57) and by the inequalities (6.53) to (6.55). 

We assume single phase fluid flow to be stable satisfying 

the condition (6.58) and particle mass density to be smaller 

1 ¥ 

than that of gas density so that - - > 0 . We are therefore 

interpreting the results for light dust with a base flow which 
in the absence of dust is stable. The shaded regions in jfig .i] 
(~ , K--) plane indicate stable regions . hotted lines, 
broken lines and solid lines bound regions which are 
respectively provided by handapat and Gupta and by the present 
author Q‘(6.47), (6 .53) to (6.55)]]. 

It is clear that sufficiency region obtained by the 
present author is much enlarged region. 
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Since for the light dust, we have not shown stable base 
flow to be always stable , the present analysis does hint 
that light dust is destabilising in nature , Hov/ever, it can 
not be regarded as a conclusive statement as theorems 
obtained are only sufficiency theorems whose violation does 
not necessarily mean instability. 

6 .8 Conclusions 

For a constant number density two phase flow, we are 
able to delineate stability regions for light dust with 
stable base flow. This region is much enlarged as compared 
to the one obtained by Demdapat and Gupta. 

A uniqueness theorem for certain initial boundary value 
problem for a constant as well as variable number density case 
has been established. It has also been established that the 
present analysis will not yield a sufficiency theorem for 
stability for variable number density case of two phase flow. 




FIG. 1 STABILITY REGIONS 
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